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ABSTRACT 

Pressure Distributions t Added-Mass, and Damping Coefficients 
for Cylinders Oscillating in a Free Surface 

A linearized theory is developed for the pressure distribu- 
tions, added-mass and damping coefficients for horizontal cylinders 
oscillating vertically with small amplitude while semi-immersed in 
the free surface of a fluid of uniform depth. The results of sam- 
ple calculations made by digital computer are presented. The cal- 
culations of the pressure distribution and the total vertical 
force for the case of a circular cylinder in fluid of infinite 
depth are compared with values measured in an experiment. 

The problem is formulated as a linearized boundary-value 
problem in the theory of infinitesimal surface waves. A velocity 
potential is synthesized from appropriate functions. The boundary 
conditions on the free surface and on the bottom and that the sur- 
face takes the form of outgoing waves far from the body are satis- 
fied exactly. The boundary condition on the surface of the cylin- 
der is satisfied at the rest position of the cylinder by an expan- 
sion in non-orthogonal functions. Convergence in the general case 
is assumed. 

The expansion coefficients for the case of a circular cylin- 
der in fluid of infinite depth were calculated by two different nu- 
merical procedures. One of these methods was chosen for further 
development to provide sample calculations for elliptic and more 
general cylinders of approximately ship-like cross-section. 

i 

The pressure fluctuation at several locations on the surface 
of a circular cylinder and the total vertical force required to sus- 
tain vertical oscillations were measured. Two-dimensional conditions 
were simulated. The measured values and values based on the sample 
calculations show similar distinctive behavior. 

The results of a general study of instrumentation for the 
measurement of pressure fluctuations in similar experiments are pre- 
sented. 
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x o Introduction 

The study of the hydrodynamics of cylinders oscillating in a 
free surface is of interest for several reasons and has therefore 
received attention for some years ® Ursell Cl9^9J formulated the 
problem of the horizontal circular cylinder, semi-submerged and 
undergoing forced periodic, heatring oscillations 9 as a boundary-value 
problem in the theory of infinitesimal gravity waves. He completed 
the rigorous solution in principle and presented the results of a 
numerical calculation that leads to added-mass and damping coeffi- 
cients® One might say the modern history of the problem begins with 
this solution by Ursell® 

Grim Cl953] follows a somewhat different procedure in synthe- 
sizing a potential solution but with a restricted number of coeffi- 
cients so that the results are approximate 0 He compares his circular- 
cylinder results with Ursell 9 s calculation® Grim extended his method 
to present approximate results for certain elliptic cylinders and for 
certain cylinders of somewhat ship-like cross-section called "Lewis 
Forms 95 in the field of naval architecture® 

Tasai [1959] returned to the method of Ursell f using in prin- 
ciple an unlimited number of expansion coef ficients, and, by a mapping 
procedure 9 extended the results of Ursell to include elliptic and 
Lewis-form cylinders® He presents the results of numerical calcula- 
tions for the added-mass coefficient and waveheight ratio (equivalent 
to damping factor) for certain elliptic and Lewis-form cylinders® 

Tasai compares his results for an elliptic cylinder (horizontal- to- 
vertical-axes ratio of 1*5) to the results of Grim® Grim’s approximate 

* References identified in the text by author and date are collected 
at the end® 
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results are markedly low for the added-mass coefficient * There is 
closer agreement for waveheight ratio* 

The work of Ursell, Grim, and Tasai cited above is restricted 
to water of infinite depth* Yu [i 960 ] added the condition of uniform 
finite-depth to the problem as formulated by Ursell • The potential 
functions used by Yu are those appropriate to infinitesimal gravity- 
waves in water of finite depth* Yu completes the solution in princi- 
ple following the method of Ursell for a circular cylinder* He pre- 
sents the results of a numerical calculation for the added-mass coeffi- 
cient and for the waveheight ratio* His result recovers the solution 
of Ursell as the water depth approaches infinity, practically speak- 
ing, the numerical results are essentially the same when the water 
depth is equal to ten radius-lengths* 

Yu also reported an experiment with a semi-submerged cylinder 
in simulated two-dimensional conditions* He measured the waveheight 
ratio and shows substantial agreement with his calculated results* 

The few other reported experiments simulating two-dimensional condi- 
tions are not readily compared* Holstein [1936] used a prismatic 
cylinder of somewhat sh5.p-like cross-section, nearly rectangular* 
Schuler [1936] used a similar form but in free oscillation, a very 
different hydrodynamic problem* Dimpker [ 193^0 used several pris- 
matic forms, including a circular cylinder, but again reports mea- 
surements for free oscillations* 

There exists work on the similar problem in three dimensions, 
on approximate methods for solving the three-dimensional problem by 
drawing on two-dimensional results, and experimental work with models 
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not simulating two-dime.nsior conditions. Including the third dimen- 
sion is a sufficient departure from the objective of this work that 
pursuit of that worthy problem is not undertaken here. 

The present paper considers the case cf the more general cylin- 
ders in a fluid that may be of uniform finite depth. The problem is 
formulated as a boundary-value problem in the theory of infinitesimal 
gravity waves, Within this framework the problem is linear and the 
solution is obtained by superposition of appropriate velocity poten- 
tialso Following Ursell, the solution is obtained in principle employ- 
ing an unlimited number of expansion coef ficients 0 For the special 
case of the circular cylinder, one recovers the solution of Yu. For 
the special case of infinite water depth, a solution is obtained which 
includes that of Tasai and is also explicit for more general cylinders 
than he treated (for example , the three-parameter cross-sec tion family 
of Landweber and Maeagno [193910 o Finally, for the special case of 
infinite depth and a circular cylinder, the parent solution of Ursell 
is recovered® 

The hydrodynamic pressure distribution for the cases of the 
previous solutions can oe derived from the wo rk of the authors cited® 
However, none state this result or present sample calculations for 
their cases. The pressure distribution is considered here in order 
to provide possible insight into similar problems of naval architectur- 
al interest and to provide calculations for comparison with experiment- 
al measurements. 

The results of numerical calculations for the case of a circular 
cylinder in fluid of infinite depth are presented® The numerical 
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procedure differs from that of Ursell for the same case. The force 
and waveheight ratio results are in agreement . The pressure distribu- 
tion was not previously reported® 

In order to compare the theoretical results with experimental 
values, measurements were made on a circular cylinder in simulated 
two-dimensional conditions® Total vertical force and water pressure 
fluctuation on the model surfacee were measured. The measured values 
compare favorably with the predicted values. Over a small part of the 
frequency range used, a distinct departure from two— dimensional condi- 
tions was observed and was due to the particular experimental arrange- 
ment® In this small range, measured values consistently deviate from 
predicted two-dimensional values® 

The results of numerical calculations for elliptic cylinders 
and certain cylinders of ship-like cross-section in fluid of infinite 
depth are presented. Certain added-mass and waveheight ratios of this 
group can be compared with those of Tasai 9 who used a different numeri- 
cal procedure. The pressure distributions and the results for certain 
cross-sections have not been previously reported® 

Numerical calculations following this procedure for circular, 
elliptic and more general cylinders in fluid of finite depth are to be 
done at the Computer Center^ University of California® The calculations 
reported here were done at this facility. 
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lie Formulation of the Problem 



2*1 Description* 

We consider cylinders of a certain class to be defined more 
precisely in the next section. The cylinder is immersed in a pre- 
viously-undisturbed fluid which may be of uniform finite or infinite 
depth. The axis of the cylinder is in the plane of the still free 
surface. Consider now that the cylinder is forced in vertical sim- 
ple harmonic motion and that steady-state conditions are attained 0 
We assume two-dimensional conditions, and that the amplitude of 
oscillation is small compared with the cylinder beam. In the steady 
state, the forced vertical motion generates a surface disturbance 
that at a distance from the cylinder takes the form of uniform pro- 
gressive gravity waves whose amplitude is proportional to the oscil- 
lation amplitude. 

We require the pressure distribution on the surface of the 
cylinder, the vertical force required to sustain the oscillations, 
and the ratio of the waveheight far from the body to the oscillation 
amplitude. 

The total vertical force in the steady-state is periodic and 
includes two components of particular Interest. Because of steady- 
state periodic conditions it is possible to resolve the vertical 
force into components in phase with the motion acceleration and 
velocity. These components do not depend on the previous history 
of motion, given steady-state. It is possible therefore to define 
added-mass and damping coefficients as ratios of the appropriate 
force component to the mass of the displaced fluid that are functions 



of body geometry and oscillation frequency only. Expressions for 
the added-mass and damping coefficients are required. 



2.2 



The Cylinders. 



The cylinders admitted are those mapped by conformal trans- 
formation of circles by 

© r\ 



where* £ = i yoe 



N 

-= C+ Z a_ . 

a ^ . 2n+l 

n=0 



-iO 



-(2n+l) 




9 



z = x + ly 



ire 



-ip 



© 




The elementary case ^zn+l = 0 
provides a circular cylinder* The case 
N ~ 0 yields elliptic cylinders 0 Lewis • 

[1929] introduced the use of this transform 

with N = 1 for a two-parameter family with a^ f a^ selected to 
provide cylinders of more or less ship-like cross-section* The 
solution of Tasai [19593 employs this two-parameter family. Landweber 
and Macagno [19593 use this transform with N = 2 and select a^, a^ 
to provide a more general selection of ship-like cross-sections* 
Prohaska [19^73 treated sections derived by transformations employ- 
ing a^ 9 and a^ f a^ and a^ t a^. 

In this work we leave N unspecified in the solution of the 
problem, and select various a^ for sample numerical calculations 

* The and y-axes are chosen positive down from the free surface 

and the angular coordinates 9, j3 are measured from the vertical 
to facilitate reference to related work cited. 
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2. 3 Infinitesimal Surface Wave Theory. 

The problem is solved within the framework of infinitesimal 
surface wave theory. We may apply this theory if we assume that 
the fluid is incompressible, that the flow is irrotational , and 
that the motion amplitudes and velocities are all sufficiently small 
to reject all but linear terms. In this problem we also neglect sur- 
face tension. 

The process of linearization is specifically employed three 
times: First, in linearization of the free -surface boundary condi- 

tion; second, in satisfying the boundary condition on the surface of 
the cylinder at the rest-position; third, in the approximation of 
the hydrodynamic pressure by using the linearized form of Bernoulli’s 
equation. These are well-established elements of the general theory 
within which we work. For a more complete discussion see Wehausen 
[i960] or Stoker [195?]. 

2.4 The Boundary -Yalue Problem in Potential Theory, 

Given the assumptions above 9 we can state the boundary-value 

problem in potential theory® We seek a velocity potential <p which 

2 

is a solution of the Laplace equation V <p - 0 in the fluid domain 
and which satisfies the following boundary conditions® 

(i) The linearised £re@-surface condition: 

K<P 4 . |2. = 0 

sy 

on the surface y = 0 outside the body of the cylinder. Here 
K = 0) 2 /g. 
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(ii) The normal velocity at the lower boundary of the fluid 
vanishes : 




on the bottom y = h. 

(iii) The radiation condition: The disturbed surface takes 

the form of regular progressive outgoing gravity waves at large 
distance from the cylinder. 

(iv) The normal velocity of the fluid at the surface of 
the cylinder is equal to the component of the cylinder’s forced 
velocity in that direction. 
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III# So?.ut±on of the Problem 



3#1 Description. 

The problem has been stated in Chapter II 0 The solution is 
described in this section and outlined in this chapter. The com- 
plete velocity potential 9 (and the corresponding stream function 
y) describing the fluid motion will be synthesized as the sum of 
appropriate potential functions. Each component function individu- 
ally satisfies the Laplace equation and the boundary conditions on 
the free surface and on the bottom. The asymptotic behavior of the 
component functions for large distance from the cylinder correctly 
contributes to a potential appropriate to outgoing surface waves, 
and thus the radiation condition is satisfied. In order to satisfy 
the last boundary condition, concerning the normal velocity of the 
fluid at the cylinder boundary, an expansion in non-orthogonal 
functions is required. It is shown that an unlimited number of 
coefficients can, in principle, be determined. Convergence of the 
expansion, proved for the circular cylinder by Urseli [ 19 ^ 9 * 19533 9 
is assumed for the general case. Thu® the complete potential and 
stream function are known and the problem is solved. Having the 
potential, we determine the pressure distribution and, by integra- 
tion, the vertical force required to sustain the oscillation. The 
vertical force is resolved into components in phase with the ac- 
celeration and velocity to give added-mass and damping coefficients. 
The ratio of waveheignt to oscillation amplitude is available from 
the asymptotic behavior (progressive surface waves) of the potential 
and from the damping coefficient 0 
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3o 2 Notation of the Geometry* 

The s- or physical-plane geometry is conformally mapped from 
the reference £ -plane by the transform 



2L 

a 



N 

C + I 

n~0 



a 2r«-l 



•(2n-fl) 



where a and a 2 n o are real so that the z-plane coordinates mapped 
from £ = ±pz“^ are 



. . =18 

z = x + iy = ire , 

N 



x = ajjosin © + £(=l) n — sin(2n+l)©J 



n=0 p 

.n+1 “2n+l 



y = afjocos © + Z(~l) n+ ^ ccs(2n+l)oj . 

n=0 p 



Special notation will be convenient for coordinate values on the 
cylinder surface- mapped from the reference circle p = p Q . 

x( body ) = aX , X = X(p o , ©) 

x(body,9 = Tt/2) = aG , G = X(p Q , n/2) 

y(body) = aY , Y = Y(yO o , 9) 

The product aG is the half-beam, b, of the cylinder at 
the free surface. The fluid depth is uniform at y = h. 



3.3 The Component Potential Functions. 

In the following sections we shall develop the complete veloci- 
ty potential <p and stream function y in the form 
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< 4 > £ cz 

B <T- 

sin ct -6* cos cot «$* p 



Y. 



a 






w 

&x& , 2 m , . * 

. os wt -f sxn oof \ + 

\ w 1 



2m 



2m 



cz 



/ ^ 



^2m( 



2m 



m~x 



\ Y 2m 



w 

4 2m 

>in oc-t « w cos cot 



Y 



2m 



Components are named as follows 



q> g source potential! 



ip sfanding«wave potential^ 

c 

^2m s ^2m expansion coefficients! 

q> 2 m multipole potential! 



W 

associated standing~wave potential, 



' <D 
' S 



The combined components ( sin cut + ' cos <*>tl and the combined 

\ T a Y*. / 



/ ^2m Ces 0)4 #2m sin ^ \ 

components ( . . ±, ^ . ) display asymptotic behavior 

\ Y~ sm o>t Yo cos out / 



appropriate to outgoing surface waves. The multipole 



<P 



2m 

r 2m 



in turn 



consists of two parts, each with several terms. 



3.31 The Multipole Potential <p ^ . 

The mu.Xti.peXe potential is written in two parts 

- ®2nT“ 9 + < p2m^ x '’ y9ll ^° ffi“l 9 2 9 3 9 ... 

The part <f‘ ’r,B satisfies the free~surfaee condition and 
2m 

the bottom condition for infinite depth but does not satisfy the 
bottom condition for finite depth. The part <p" m (x,y ,h) is added 
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to make the total satisfy the finite-depth boundary condition, follow- 
ing Thorne [19533* The function cp 2 Q (x,y,h) independently satisfies 
the free-surface condition and for h oc, ^ 2 m ' x *y «**) -*• 0 * 

The index 2m will be understood when it is dropped from the 
symbols and . 

(a) Part 



,, a^ m cos 2m|3 f , s n a^ m+ ^ n+ ^^cos(2m+2n+2) p 

<f(r ' 6) ' E+ I 7 /- 1) a 2n+l Z 

r n=0 r 



K 



2m 



2m- 



a^cosUm-Dfl ^ K f M) » 

1 - 2 "'- 1 fed 



a- , 2m+2n+2 

2n+l a cos(£m+2n+l) p 

2m+2n+l 2m+2n+l 



To show that tp‘(r,0) satisfies the free-surface condition, 
it is convenient to express that boundary condition in another form: 



^ * |*' + Kcp” + |*b = 0 , p = 9 = n/2 . 



2m dy 



dy 



3y 



The part K(p' + 5<p'/dy , independently zero, is written 



K<p‘ " r If = °* P = n/Z ' 

When the substitution is made and evaluated at (3 = n/2. f the remain- 
ing terms are 



a^ m [cos m n + sin( 2m-l) 7i/2]/r^ m + ^ ( -l) n 

n=0 



2s r 

a^. . a Lcos s n + 

dn + 1 



+ sin(2s-l)n;/2]/r • 



We have introduced 2s = 2m +2n+2. Now, using the identities 

cos mn = (-l) m and sin( 2m-l )n/2 = (-l) m+ \ it is clear this sum 

is identically zero; therefore, cp 1 (r,[3) satisfies the free-surface 



condition. 




* 
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The function tp^ (r,{J) satisfies the bottom condition only 

in the case of infinite depth. The total cp_ = cp' + cp" will be 

dm dm dm 

adjusted to satisfy the bottom condition at finite depth in (c) 
below. 



(b) Part <p 2 m ( x »y«k) 



CD 

tp"(x,y,h) — / [c_(k)sinh ky + c^(k)cosh k(h-y)]cos kx dk , 
o 



where c^(k), c^(k) are found below and assure convergence. 

Carrying out the operation (K + d/dy)cp'' leads to the 
definition 



c 4 (fc) = - Cj(k) 



kh - ksinh kh 



to assure that tp"(x,y,h) satisfies the free-surface condition 
at y = 0 . 



(c) The Finite-Depth Boundary Condition for the Multipole. 

The coefficient c^(k) in the integral term tp"(x,y,h) 
is determined by using the bottom condition, 



dtp 



2m dtp* dtp n _ . 

, — = “ = 0 , y = h 

9y 3y Sy 



To find the first term, use the expansion (for y > 0) 



cos 



n 



/ k n **^e“^coskx dk . 
• o 



[Whittaker and Watson, A Course of Modern Analysis, Cambridge 1927. 
para. 12.2 Ex.2], With this expansion, one may write for tp',when 



y > 0 , 



14 



GD 



"ky 



V* = / L(k)s“ ‘ coskx dk 
o 

N 



where 



L(k) = 



, 2m-2 2 m,.... , ^ 
K a (k-nf.) 

(2m-l) i 



2s«2 2s, 



21-0 

The bottom condition is satisfied when 



■'T ,■ - \ k a ( K+ k ) 

f „A“ 1; 



3<p 



2m 



QD 



y»fc 



-kLs ' + C-,(k.) k cosh khjcoskx dk 



■kh 



vanishes. The definition c^(k) = Ls“ "/cosh kh assures this. 
With c^(k) determined, cp”(x,y,h) is written as a Cauchy 
Principal Value integral (/)„ 



, , N , , -kh K sxnh ky - k cosh kv 

<p_(x,y , h) = / Le . “ — — coskx 

5 ,J ’ K cosh kh - k sxnh kh 



dk 



= / 



oo -kh,... . v, 2m-2 2m „ . . , , ,, 

s iK+k;k a K sxnh ky - k cosh ky 

(2m-l)J K cosh kh - k cosh kh COS * x + 



, CQ s'° K:!l (K-t-k)k ? ~'~‘~a^ ~ K sinh ky - k cosh ky 



3»32 The Combined Mult5.pole and Associated Standing-Wave Potentials. 

The radiation condition [2.4{iii)j requires that the asymptotic 

behavior of the potential for large x and y = 0 be appropriate to 

outwardly progressing surface waves. The part <p' of the multipole 

tends to aero for j x j However, the part <p ? of the multipole 

does not tend to zero for ixj y = 0. The associated standing- 

w 

wave potential tp^ is shown below to combine with the multipole to 
form a total which satisfies the radiation condition. 

Since <p'(r,£) 0 for jx f “> ®, cpp m (x,o,h) ^'^"(XjOjh) for 

| x | ^ a « It is shown in Appendix A-l that the following asymptotic 



expansion holds; 
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cosh K (h-y) 

II , , N _ o J 

<p_(x,y,h) ~ - E 



0 „ , , sin K x as 

2m cosh K h oil 



|x -> CD , 



and 



(pp(x,o,h) ~ - E 2m sinKjx 



as x -» oo , 



where 



E- (K a,K h) = 



2ti 



2m o ’ o 2K h + sinh 2K h 
o o 



» 2m 2m 
K q a 

(2m-l) * + 



N 



Z <-i> 



n=0 



„ 2s _ 

a_ .. K 2s 
n 2n+l o a 

(2s-l)» 



and K Q (K,h) is defined as the positive real root of 



K = K tanh K h . 
o o 

For finite h > 0 the asymptotic behavior of (p"(x,o,h) 
as | x | ->■ gd with an appropriate harmonic function of time, is 
identified as typical of a standing wave in a fluid of depth h. 
We add a standing-wave potential consisting of the function 



w 

<P 2m = c 2m cosh K o (h-y) cos K q x m=l,2,3,... 

multiplied by a harmonic function of time. The coefficients 
are chosen in the following so that the combined potential repre- 
sents progressive waves and the radiation condition is satisfied. 

To determine c_ , consider 
2m' 



*■ 1 



. art 



<P 2m cos cot + cp^sin cot] ^ cp 2m (x,o,h)cos cot + <p 2m (x, o,h)sin cot as |x| ->oo. 



y=o 



The right-hand side is asymptotically equal to 



- E 2m sin K o |x| cos cot + c ?m cosh K^h cos K n x sin cot as |x| ->> 00 . 



"2m 
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By defining 



c.» ~ E~ ^osh K h 

2m ills' o 



the asymptotic behavior of the combined functions is 



- E 2m Sln (K o H “ Wt) 9 

which displays the required behavior*. 

We also form the sum with conjugate harmonic time functions, 



which for 



*2m 

CD , 



W 

sin cot - <f> 2 m cos d)t , 
y = 0, tends asymptotically to 



- E 2m cos (K^jxj - cot) „ 



A linear combination of multipole potentials, adjusted as 
above to be progressive°wave producing, is 



CD CD 

Z p 2m ^2m Cosajt+ '4n sina)t)+ Z q 2m ((p 2m sin ^ “ ^2m C0Sa)t)o 

m=l m-1 

The coefficients p_ (Ka) and a_ (Ka) are to be determined, and 

^2m *2m * 

are assumed to behave so that convergence is assured,, 

The asymptotic form of the free surface y = 0 for large 
x > 0 is 



1 d W w 

t| = “ ~ “ 7 P-, E cos (K x-tot) = “ ^ q, E sin (K x-oot) „ 

g dt g *-.^2m o g jf- *2m o 

In complex notation, if we reserve the operator i for 
geometric variables and introduce the operator j for functions 
of time t, this sum is 
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CD 



Ee d 2 Cp 2m + 3 * 3 



W . ~jO)t 






We will later make use of the corresponding stream functions 



w 



Y 2 m and y 2 m * which may be found from 



f 2m (z) ” *2m * ±Y 2m 



N 



= (-I)". 2 "*- 2 " . Z !-l)a (-l)S- 2s -(-l)” ji£ r a 2 ”z- (2n - 1) 

n-0 + m “ 



^ / , N n / , N s iK 2s -(2s=l) 

2 /- 1 ) a 2nal ( - :L> 25ZT a z 

n=0 



and f™ (z) 
2m 



? L £ ~ kh [iK sin kz ^ k cos kz 
K ccshkh- k sinh kh 



w . w 

him * 1 Y 2m 



- c_ cosh K (h+i z) . 
dm o 

With these definitions* the complete wave-producing multipole is 



GD 



Re, I (Pp m ♦ ^ W (f 2m Ca) + j f 2m (z))e “ jWt 






3.33 The Source Potential. 

i 

The potential function that represents a source at the origin 
and that satisfies the free-surface boundary condition is 

cd -'ey . 7 , qd -mx , v 

, e "coskx ,, =Ky . Tr e (mcosmy=K sin my j 

£- - - dk = its J sxnKX" / — — ^ ^ ^dm . 

o o m + K 



It is easily seen from the form on the right that this satisfies 
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the f ree-surfac e condition. The bottom condition, however, is not 
satisfied in the case cf fluid of finite depth. To satisfy the 
finite-depth bottom condition, we add a term in the manner of 
Thorne [1953] , 

/ [c^ (k)sinh ky -u* C 2 (k)cosh k(h-y) ] cos kx dk . 
o 

The coefficients c^(k), C 2 (k) are determined below and assure 
convergence . 

This integral term satisfies the free-surface condition 
given that 

C 2 (k) - - c^(k) ^ cosh kh - k sinh kh 

The coefficient c, (k) is determined so that the combined 

term 



^coskx j, ® N r • , , K coshk(h-.y) . , 

£ =— : dk -v £ c, (k) Lsinh ky - . , . ' . , v ■ ; — r-. - , cos kx dk 

* K-k 1 J K cosh kh-k sxnh kh 

o o 

satisfies the finite-depth boundary condition. This boundary 
condition, d<p/dy = 0, y = h, leads to 

-kh ~ 

c (Jr) - — „ 

l'- ; K-k cosh kh 



With this coefficient, the source potential is 



tp 



. coshk(k°y) coskx .. 

^ K cosh kh-k sinh kh 
o 



When h -¥qb, this simplifies to the source potential 



given first 
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3o3^ The Combined Source and Standing-Wave Pctentialso 

The radiation condition is satisfied by adding a standing- 
wave; potential to the source potential 0 with the necessary harmonic 
functions of time 9 to form the term 



<p cos (ot 4* <s> sin out 

© ■ s 



The standing-wave potential function is 



cp = a cosh K (h-y) ©os K x 9 

©CO 0 



with a Q defined as given below* 

Appendix A-2 shows that the asymptotic behavior of the source- 



function <J> for jxj a? 5 



0 is 



E sin E ix 1 

s o 1 ! 



where 



2ix( cosh K h) 

E fK — P 

s* o ' 2K h. + sinh 2K h 



Now a is defined as 
o 



E 



a = 



o cosh K h 

© 



so that the asymptotic fcrm 9 as ix: =# ao , 

(<p cosa)t+ cp sinajt) ** e. cosh K h cos K I xj cos ut + 
c s \ - c o © 1 • 

. 17 “ o 



+ E sin K he sin cot 

& o 1 1 



may be written 



E coe 'K |xi - a>t) 

SI Cn ' 
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The source potential and its associated standing-wave potential 

0° coshk(h-y) cos kx ,, ^s , ,, ,, > 

/ t-TT 1 ■ — r rr dk sin cut + r - r cosh K (h-y ) cos K x cos U)t 

K coshkh - k sinh kh cosh K h o J o 

o o 

therefore satisfy the free-surf ace , finite-depth, and radiation bound- 
ary conditions. 

In the special case h -> ao this may be written 

r co -mx, ir . 

-Kv . „ , e (m cos my - K sin my , \ ^ -Ky 

7te ^sm Kx = / sj-“ p ‘‘-dm sin cut + ne J cos Kx cos cot , 

L o m + K J 

which is the form used by Ursell [19^93 and Tasai [19593. 

The form of the free surface for large x > 0 for either 



case is 



n , CUE 

ldtp s . f Ir 

ri = — rr ~ sm ( K x - cut ) . 

' g 3t g o 



For the special case h CD , E g it 



The stream function corresponding to this potential may be 
found from 



Re ( f ( z ) + jf (z))e' 

J c s 



■jtut 



where 



f c ( z) = cp c + iy c 



2 n cosh K h 

rr— ; : — : — rrr. r cosh K (h+iz) 

2K h + sinh 2K h o' 



f s U) = * s + iY s 



oo cosh K (h+iz) 

$ 2 



dk 



K coshkh - k sinh kh 
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3.4 The Potential Solution. 

The potential solution written as a function of the complex 
variable z = x + iy is 



The coefficients P 2 ffl (Ka), are determined later by use of 

the normal-velocity boundary condition. The right-hand side has 
been normalized by the coefficient J(K Q h) in order to display the 
surface waveheight h^ at large distance from the cylinder. 

The function J(K h) may be determined as follows, for exam- 

o 

pie, from the asymptotic development of the potential for large 
x > 0, y = 0. 

gh w p oo 

cp(x 9 o,h) ~ — — E cos (K x-Cjot) + 2L P 3 E. sin (K x-ajt) - 
J %0)J L s o _ 2m 2m o 



The corresponding asymptotic surface waveheight as x cp is 








This may be written 



sin (K c x - wt + x) 
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when we define 




2 




For h ql> , we recall that E.^ ■*> C , E^ it , so that J ^1 , 
All quantities in the potential and stream function 



except the expansion coefficients for a given cylinder, water- 
depth, and heaving oscillation have been determined,, When the 
expansion coefficients are determined in the next section it will 
be convenient to have the .following notation,. The subscript a 
designates evaluation on the cylinder boundary. 



Y 




L Y g 





-?] 



body 





a- 

‘dm ' 
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y] - —— - (D s^n cot + C coscot) 
, , tcox; 

body 



CD CD 

W 



D(Ka f 0) = y + 7 P, + 7 q, y, 

9 sa 2m * 2ma ^-_ H 2m l 2ma 



m=l 



m=l 



CD 



CD 



w 



C ( Ka 1 9 ) Y ea + Z p 2m Y 2 ” Z q -2m *2raa 

m=l ma m=l 



In addition we designate the evaluation of y] for 9 = 0 = n/2 by 

body 

gh 



yJ 



w 

TtCOJ 



body 9 it/2 
B(Ka) = D ( Ka , tc/ 2 ) = y 



( B sin cot + A c o s cot ) 



CO 



y W y 

sa n/Z ^-. P 2m ^2ma n/Z * q 2m Y 2ma n/Z 
m=JL m=l 



QD 



A(Ka) = C(Ka,n/2) - Y ca ♦ Z P 2m V 2ina „ /2 - I % r” na „ /2 . 

m~X m=l 



The apparent redundancy of notation has the purpose of facilitating 
reference to related solutions previously mentioned. 



3o5 Evaluating the Expansion Coefficients. 

The vertical velocity component that is perpendicular to 
the cylinder surface must equal the normal velocity of the fluid 
at that point. This boundary condition is to be satisfied at the 
rest position of the cylinder by suitable choice of expansion co- 
efficients p_ (Ka). q_ (Ka). 

2m y 2m 

The vertical motion is y = h, cos (cot + e). After using 



the substitutions 
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dx _ d%_ 
ds ~ dn 



cos a 



- . ix-i 

dn ds J 



body 



the boundary condition is 



_ dy _ dy dx -, 
ds - dt ds J 



body 




y v 



After integrating and setting the integration constant equal to zero, 



«-£*] 



dt 



body 



We evaluate this expression at the special location p = Q = it/2,, 
using the notation of the preceding section: 



y] 

body , rt/2 
Then we may use 



dt ~ nqjaGJ 



gh 



w 



TtCOJ 



(A cos cot + B sin art) = 






(A cos (ut + B sin art) 



to evaluate the same expression in general locations on the cylinder; 

SK 



y3 



body 



(C cosut + D sin art ) = - x] 



body 



gh, 



w 



(A cos art + B sinart)aX. 



itoraGJ 

Equating the coefficients of sin art and of cos art, we obtain 



C(Ka ,9) - 



X(Q) 



A ( Ka ) = 0 , 



D(Ka,9) - Up- 



B(Ka) = 0 
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Substitution of the definitions of A, B, C, and D and introduc- 
tion, for convenience, of the notation 

f 2m ( Ka ,e) = - r 2man / 2 - r 2ma » 

e ( Ka Q) = - v W - v W 

g 2nr ’ G ^2ma n/2 ^2ma ’ 

leads to the two equations 



CD 


CD 


X 

G ^can/2 ’ 


£ P 2m f 2m " 
m=l 


^2m g 2m ~ ^ca 

m=l 


00 


CD 


v 


^ P 2m g 2m + 
m=l 


^2m^2m ~ ^sa 

m=l 


A 

G ^san/2 



These two equations display the roles of p 2m (ka) and q 2m (Ka) 
as expansion coefficients. In a given example the coefficients 
can in principle be evaluated. In practical calculations some 
sufficient finite number of coefficients can be found by the ap- 
proximation methods employed for expansions in non-orthogonal 
functions, as in the sample calculations given later. 

In the special case of fluid of infinite depth, g 2m (Ka,©) 
is identically zero. 



m— 1 ,2,5,..., 



m— 1 ,2,3,..., 



3.6 The Ratio of Waveheight to Oscillation Amplitude. 

We return to the expression for the normal velocity condi- 
tion evaluated at (3 = © = m/2 given above. 



gh 
B w 

nooJ 



(A cos cot + B sin cot) = 



" ft (aG) 



= coh b sin (cot + e) . 
o 
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The trigonometric expansion for sin (out + e) leads to 



A 



B 



2, h 

0) b o T . 

- — jt J sin e 

g h 
° w 



2 k h 

0) b o T 
— — ~ n J cos s 
g h 
° w 



These give, respectively, by division and addition after squaring, 

, =1 A 

e = tan - , 

h w _ Kb k J 

h o ft? + 



The first result defines the previously -unspecified phase angle 
e of the motion. The second gives the ratio of the height at 
infinity of the generated progressive waves to the oscillation 
amplitude, 

3,7 Pressure on the Cylinder Surface, 

Gauge pressure at a point on the cylinder surface with 
reference to the constant pressure on the free surface is 

(P - P Q ) - pSJ - p |f 

from the linearized form of Bernoulli’s equation. The first term 
on the right contributes a hydrostatic pressure due to the loca- 
tion of the point at rest, 



7 * aY , 

and a fluctuation of the hydrostatic pressure due to vertical motion 
y 0 gh Q cos (art + e) , 



The hydrodynamic pressure is 



p ( Ka , © , t ) 



- F 

/Dgh 






w 



teF 



[ - M COS Cu't + 



N sin cot ] 



where we have used the notation for the potential tp introduced 
in Section and M(Ka s ©) { N(Ka,G) are to be evaluated at the 
point on the cylinder of interest., 

The acceleration of the given motion 



y = h cos (cot -fs) 
o 



h. 

/a 2 + B 2 



(B cos art - A sinoot) 



is 



t: 2 h 



S - 



VA 2 -i- 3 2 



(A sin a:t - B cosoot) 



The hydrodynamic pressure p is resolved into components in phase 
with the acceleration and with the vertical velocity s 

P g '\ KB - 5 * NA 



P ^ Ka o G » ^ ^ f -> 

A + 3 



(A sin cot ~ 3 ccs cot) + 



jP gn w MA N B 

vF 2 2 

" A + 3 



(~A cos cot •* B sinoot) 



After substituting for the wave-height h^ , this may be written 



MB + NA 



p F b ~~~J ? + P h ~ 2 



A” •> B 



MA - NB • 

coy 

A“ + 



or in units of the hydrostatic fluctuation yOgh^. 



9 
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pk h 



- Kb 



MB NA 



y 



+ Kb 



A + B ' <S h 



MA 

,2 



NB 



B' 



2 o,b 






where M(Ka,©), N(Ka 9 0) are evaluated for the location of interest. 
The total change in fluid pressure must include the fluctuation in 
hydrostatic pressure due to vertical motion. This may be added in 
normallized form y/h. to the above expression to give the total. 
Examples for specific cylinders are given later. 

5.8 The Vertical Force Required to Sustain Oscillations. 

One parr of the vertical force is that opposing the result- 
ant of the hydrodynamic pressure. This is called the hydrodynamic 
force. The hydrodynamic force per unit cylinder length is found 
by integrating the vertical component of the hydrodynamic pressure 
acting on an element of area. 



dF = p(Ka,0,t) cos ads . 



After substituting cos a = dx/ds and also 



dx 

ds 



ds = dx = a 



cos 0 + 



N 

z< 

n~0 



■X) J 



(2n+l)a 



2n+l 



A 



2n+l 



cos (2n+l)Oj dQ 



- aW ( 0 ) d© 



in the form 



v W (Q) 

dx - b d© 

U 



where b - aG is the half-beam at the free surface 9 the force 



becomes 
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x/2 , „ 

7 = 2b f p ( Ka , © , t ) d© „ 

U 

O 

Recalling the notation for the pressure introduced in the last sec* 
tion 9 one may express this in the form 

2 /' gh w b 

F = - — —r — - (N sin <ot - M cos cot) 

71 J O O 



where 



u/2 



W(S) 



M o (Ka) * / M(Ka,0) d© , 



»i/2 



N (Ka) 
o 






N(Ka,©) d o 



The hydrodynamic force is resolved following the procedure 
of the preceding section, 

?.p ghb M.BeN.A 

F(Ka,t) - (a sin cot - B coscot) 

A + B 



?.pgh ’o M A~N B 

w ° o , . _ . . 

•> ™ — — s— — (“A cos cot - B sxncot) 

no a 2 + p 2 



Substituting for the wa ' eheignt h t , one iindc 



w’ 



F - 2pb 



0 M B + N A „ 

2 a o 



- M A - H 3 , 

. 2 o o * 

coy 



5' y •> «spd 

A + B“ A 2 + B 2 



It is convenient to norma lire rhe hydrodynamic force in units 
of the change of hydrostatic, force (naval architectural ’’tons- 

J 

per-inch immersion” ) » 



_ MB+NA ” M A - N B • 

™rr” ~ Kb (-f— ) -s- Kb — %- (-f~) 

2ycgbn o a 2 B 2 w 2 b a'' + B 2 
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The total force fluctuation is found by adding the normalized 

buoyancy change, y/h , and of course account must be taken of 

o 

the inertia of the mass of the body and net buoyancy at the rest 
position. 



3.9 The Added-Mass and Damping Coefficients. 

Since the fluid motion, and thus the pressure on the cylin- 
der surface and the vertical force, is periodic, it is meaningful 
to define added-mass and damping coefficients and they will not 
be functions of time, i. s„ of the past history of the motion. 



3.91 The Added Mass and Added-Mass Coefficient. 



The added mass is the ratio of the hydrodynamic force in 
phase with the acceleration to the acceleration 



n ,2 

m = rr b p — 

2 r it 



,, M B + N A 
4 o o 



2 2 
A + B 



The asymptotic value of the added mass for large Kb in the spe- 
cial case of infinitely deep fluid is given in Section 4.4. 

It is convenient in this problem to define an added-mass 
coefficient as the ratio of the added mass m to the mass of the 
fluid displaced by a circular cylinder of equal beam. 



.. M B + N A 
w. _ 4_ o o 

TC ,2 ~ 1 J .2 -2 

p A + B 



To normalize the added mass with units of the actual displaced 
fluid mass, one may use the volume of the cylinder per unit length, 
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V - 



t: b_2, 
2 G 2 



Pr 



N (2a + l) <a 2r _ +1 )-n 



2 

n=0 



C2n+1>2 

Po 



so that the ratio of the added mass to the mass of the displaced 
fluid is 



~rf , MB-:- N A 
m o_ 

/» = v *' a 2 ; b 2 ‘ 



The coefficient defined in this manner is inconvenient for the 
case of a flat plate oriented perpendicular to the motion* In 
this case V = 0$ however* the added mass is still defined and 
finite* 



3*92 The Radiated Energy and the Damping Coefficient* 

An energy balance must exist between the average rate of 
work done in one cycle of vertical motion and the average energy 
flux radiated by the outward-going surface waves* The average 
rate of work done over one period T is 

er;uh 2 /T 
o 

where c is the factor o:: proportionality for the component of 
hydrodynamic force in pha.se with the velocity* 

F =s c(dy/dt) * 
v 

The average energy flux per unit width of periodic outgoing waves 



\ h 
2 w 



p< 



(j) 



on two sides is 
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Therefore , 



- h •' 2 - 2nJ*~ 



0 = Or) 

o 0) 



2 2 
A + B 



0) 



When this is compared with the results in Section 3*8 f 



M A - N B = ti 2 J 2 /2 . 
o o 



This is useful for checking numerical calculations and simplifies 
the expression for the component of vertical force in phase with 
the motion velocity, 



^ ^2 (n 2 J 2 /2) * 

F_ = 2 pb ~ L — coy 



v 



■r 



2 2 
A + B 



The ratio 



v 

n, 2 

z° r 



2nJ 



A 2 +, B 2 



coy- 



leads to the convenient use of 



2nJ 



2 2 
A + B 



as a damping-force coefficients If the added-mass coefficient is re- 
ferred to the actual mass of fluid displaced by the volume V, the 
corresponding definition of the damping-force coefficient is 



yOV «(0$- 



n y2 

2^ 

V ft 2 2 
A + B 



2nJ 
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IV o The Special Case of Infinite Fluid-Depth 

The case of a fluid of infinite depth deserves special men- 
tion 0 An important reason is that the expansion to satisfy the 
normal-velocity boundary condition on the cylinder can be done with 
different functions* A second purpose is to summarize the simpli- 
fied definitions* Finally, with respect to the added-mass coeffi- 
cient, a correction factor accounting for free surface effects if 
readily definedo For practical calculations, this case is simpler 
and is of interest since the infinite-depth results are approached 
rapidly with increasing depth ratio* 

4*1 Special Definition of the Multipole Potential * 

The multipole potential which may be used in case h = oo 
is designated by anc * ^ e “^- nec ^ ^y 

(2n+l)a 2ri+1 cos (2m+2n+l)Q -i 

(2» t 2 nt D / > (2mt2n+i) y 

This replaces and s; '- nce ® as ^ ->ao, there need 

be no counterpart to (plj (x,y,h)o It is evident that van ^ s ^ es 

for |x| oo and thus has no role in satisfying the radiation bound- 

w 

ary condition. It is consistent that ^ - 0 for this case. It is 

shown below that <p_ satisfies the free-surface condition, 

2m 

The free-surface condition 

= 0 , © = + n/2 , | x | > b , 






cos 2m© 



2m 



P 



2m 



Ka 



cos{2m-l)Q 



N 

Z (-1) 

n=0 



can be written 



3k 



** * VS0SST ft * 0 



by noting that dx/d© - C when 0 = + n/2. With 



Si 

d© 



N 

p sin © + 7 ( -1} 

n=0 



n+ 



i (2 " 1)a 2n+ i . 






/> 



2n+l 



sin (2n+l)©J 






r~ N 

, - Z 

L n=0 



(2n+l)a 



2n+l 



p 



kr] 



, 9 = +_ ti/2 



the free-surfaee condition is 



^ - j.Un.1) * . f| - 0 , 



© = + k/ 2 , | x| > b . 



Simple substitution paralleling Section 3„31(a) demonstrates that 
<P 2m ♦ like <P 2m (r,{3), satisfies the free-surface condition for all 



The corresponding stream function is 



- _ sin 2mQ 

~ 2m 

P 



v Ka 



sin (2m°l?© 
_(2m-l)p 



N „ (2n+l)a 

_ f ,»n . 

^m-TT + L . =X; (2ra+2n+l) (2m+2n+l) 

n= 0 



2m- 1 sin ( 2m+2n+l)9 

P 
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The stream function evaluated on the cylinder is used in satisfy- 
ing the normal velocity condition,. The part y ] ^ es iS na ^ et ^ 



Y2ma * as before, and at © = it/2 , 



body 



*2ma,„/2 “ Ka - 1) 



m-1 



1 ? (2n + 1)a 2m.l 

(2m-l) p 0 2m "' 1 ’ r.l'o J2m+2ntl) ^( 2m * 2n+1 ) 



j 



4,2 The Source Potential and the Complete Potential Solution* 



The definition of the combined source and standing wave 
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potential (ip cos cot + (P sin cot) remains as in Sections 3.33, 3*34, 
c s 

where it was written for the special case as 



[ a> -mx., 

-Ky . xjr e Km cos my - K sin my , 

Tte J sxn K x - / — — 5-^ — -5 — 2 dm 

o 



2 TT 2 

m * K 



sin a)t + 



«Ky 

HE 97 COS Kx COS 0)t o 



The conjugate stream function y sin a)t + y cos a)t is 

s c 



[• 



.Ky 



rr CO ~mX, . Tr \ 1 

“Ky T , r e (m sin my + K cos my) „ i 

■us ,7 cosKx + / if ^ dm sin cot + 

o m + K J 



+ Tte J sin Kx cos cot . 



For this case the potential and stream functions equivalent to 
those given in Section 3°4 are 



<P _ h w g 
. ~ noo 



rr- 


C 3D 


^2m\ . 

] cos cot + 


/(p 

rs + 


CD 

f ^2m 


sinortl 


1> Y C 


m=l 


w 




m=l 


W J 



The expansion coefficients P 2 m (Ka), q 2 ffl (Ka) are to be evaluated 
below following the procedure of Section 3«5. 



4.3 Evaluating the Expansion Coefficients,, 

The normal -velocity boundary condition on the cylinder 
surface is used to evaluate the expansion coefficients. Since 
y(yC> o ) = - (dy/dt)x( p Q ) , the two equations 

C(Ka,©) - ^~-A(Ka) = 0 , 



D(Ka,9) - B(Ka) = 0 
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are found as before, but. now with 



oo 

A(Ka) - Y can y 2 > Z P 2m ’ 2man/2 

m=l 

By. Ka,. ^san/2 " Z, „ ^2m ^2mait/2 



m=l 

OB 



C ( Ka , rc/2 ) , 
D(Ka 9 w/2) , 



C(Ka,0) = T + 2. T 2ml , 

m=X 



CD 



D'Ka.S) = Y * I V 2ma . 

m=L 



With these definitions, the expansion equations are 



CD 



I Pan, f 2m * Y 



X(9) 



m=l 



CD 



2m 2m *ca G ? can/2 



^2m " 2m ^ sa * G V sa7t/2 

m=ij_ 



where is now given by 



x(©: 



f 2 m (Ka 9 ©) = Q ' Y 2ma7E / 2 " Y 2 ma 'Ka,©) 



This is the set of expansion functions which may be used in the 
infinite-depth case to replace the equivalent unbarred set 6 

Completion of the problem and the results are identical 
in form with the corresponding steps in the general case„ 



Asymptotic Value for the Added Mass and Definition of the 
Free-Surface Correction,, 

The added mass is the same as in Section 3o91, 

_ . M B + N A 

Tt , 2 H- C C 

m = 2 0 F * ~r-— • 



A + B 
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The asymptotic value of the aided mass for large Kb (in 
:his case of infinitely deep fluid) is 



m 






TCb 2 a 



2 G 



(fc> 4 to 



N (2n+l) (a_ , ) 2 l 

+ J — . ? n+1 



(2n+l}2 



O 



which may he found in a manner following Landweber and Macagno [1957] ° 
We recall that the origin of the result called m here is in the 

OD 

work of Lewis [1929]® The ratio of the asymptotic Yalue m to the 

05 

mass of the fluid displaced by a circular cylinder of equal beam 9 
2 

ub yo/2, is the inertia-coefficient C of Lewis [19293 or the added- 
mass coefficient of Landweber and Macagno [1957« 19593 » 

The ratio of the added mass m to its asymptotic value m^ 
defines a free-surface correction factor that is a function of the 
non-dimensional frequency Kb only 0 



% = 



a 

m 



QD 



(rtV2) 

m 

GD 



M B 
o 



N A 
o 



+ B 



o 



2 

The value of for an elliptic cylinder is Tib yo/2 e from the 

definition above 0 Therefore 9 in the special case of elliptic cylin- 
ders (including the circular cylinder) 9 the free-surface correction 
factor k^ is numerically equal to the added-mass coefficient de- 
fined as the ratio of the added mass to the mass of the fluid dis- 
placed by a circular cylinder of equal beam® 



Examples of the free-surface correction factor k^ were 
calculated for various cylinders and are presented in Part B of this 



work 
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The approach used to define k,^ could be used to define a 
correction factor for the free-surface effect on the added mass of 
cylinders in fluid of finite depth but the asymptotic value m^ 
for the case of finite depth is not yet simply expressed,, 

A first-order approximation of the added mass m for ellip- 
tic cylinders valid for Kb ■=#» 0 is given by Ursell [19^93 for the 
infinite-depth cases 

m --- “ t~ log Kb - log (1 + —■) + 0.23] « 

% u £ £ D 

Here d is the draft, y(p = 0), and b remains the half-beam, 
x(P = ti/ 2). It appears from various calculations that this is a 
useful approximation for other than elliptic cylinders, but this 



is not established as a theorem 
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PART* Bo Numerical Analysis and Results of Sample Calculations*, 

Vo The Circular Cylinder in Infinitely Deep Fluid 

The principle objective of the first sample calculations was 
to provide the calculated pressure distribution on the surface of a 
circular cylinder in fluid of infinite depth. In addition, the 
added-mass and damping coefficients were calculatedo These coeffi- 
cients were calculated by Ursell and the present calculation con- 
firms his results and extends the range of the frequency variable. 

Our method of calculating the expansion coefficients P 2 m (Ka), 
and of finding the added mass coefficients is notably differento 

The calculations were done on the IBM 704 digital computer 
at the Computer Center, University of California, Berkeley. 

5.1 Calculation of the Expansion Coefficients. 

5oll Choice of Method. 

The definitions of most of the functions required in the cal- 
culation of this case are simplified. The mapping is trivial ; all 

a~ ^ are zero, and p - 1 is convenient in that the unit circle 
2n+l ’ f o 

in the reference plane corresponds to a cylinder cf radius a. The 
angular coordinates are equal, B - ©. The specially defined 
of Chapter IV is identical with the general case Chapter HI* 

We will use 6 for the frequency parameter (here Ka) in all 
numerical calculations. 



The expansion equations of Section 4.3 become in this case 
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oo 

^ P 2 m (6) f 2m (6 > 9) = *ca (6 - 9) * sin 9 r ca (».*/2) 
m=l 

ao 

£ q 2m (6) f 2m (6 SSL ^ ^S3L 9 / ^ 9 

m=l 

where 

X(9) = sin© , Y(9) = cos © , 

X(n/2) = G = 1 , 

Y (6,9) = COS ® sin (6 sin©) , 

CSL 

Y (6,©) = e ^ cos ® [(^ + log 6 + J. ~7— cos n©) sin ( 6 sin ©) - 

so. £ _ n 1 n 

n=l 

co ft n 

( 9 + 2 s ^- n n ® ) cos ( 6' sin 9) 3 , 

n=l * n 

f 2 m (&»©) = - sin 2m 9 - (^m^i) [sin (2m-l)© - (-l) m_ ^ sin 9] . 

We have substituted a form of y (6,9) from Appendix A-3«^« A 

Scl 

method must be chosen to find P 2 ffl ( 6 ), ^ 2 m^^ from the expansion 
equations for m=l , 2 , 3 , <> <, » , their computation being terminated in 
practice at some finite m , The first method chosen was to generate 
m simultaneous linear equations each for P 2 m (m=l, 2 , . . .m) and 
q 2 m (m=l, 2 , o . „ ,m) by proceeding as follows,. Each expansion equation 
was multiplied by sin k9 d9 and cos j9 d9 with k=l,2,3,.». and 
j*=0, l,2,oo. terminating at k, (j -1) so that m = k + j . The m 
equations (each for P 2 m « q 2 m^ were integrated over the interval 
0 <_ 9 _< tc/2o This will be called the "method of integrationo" The 
m equations were solved for m values P 2 m « ^ 2 m the Gauss 
elimination procedure 0 Yu [i 960 ] used the method of integration 
but with cos r©, r=l,2,3, ° . „ ,m as the only multiplier in his 



calculations 
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The second method used for sample calculations was to gener- 
ate m simultaneous linear equations from each expansion equation 
by evaluating each at more-or-less arbitrarily chosen angles in 

the range 0 < 0^ < m/2 with i=l, 2,3, <,<,<> „m „ Ursell [1949] and 

Tasai [i960] chose to take values m=l,2,3, <, „ „ ,m with m < m so 

o o 

that m values cf p~ 9 q~ could be found from m simultaneous 
o dm dm 

equations by the method of least squares to satisfy the truncated 
expansion equations at © = ©_, in the least-squares sense, In the 
sample calculations reported here the m equations were solved di- 
rectly by the Gauss elimination process for m expansion coeffi- 
cients© This will be called the "discrete-angle" method. 

In all cases the calculated expansion coefficients were tested 
by substituting them in the left-hand side of the basic expansion 
equations (Section 5*11) and comparing this result with the calculated 
right-hand side at various values of © = ©__ © This calculation was 
printed from the computer with five significant figures. In the case 
of the calculations by the discrete-angle method, test values 9 = ©^ 
that equal the discrete-angles ©^ lead to exact agreement and gave 
proof the matrix was properly inverted© For interstitial 9^, the 
difference in right- and left-hand sides is an accuracy index not 
generally zero© This index would not 9 in general 9 be zero for any 
0 < 9, < ti/2 la the case of expansion coefficients calculated by the 
method of integration from integrated equations and substituted into 
the basic expansion equations.. 

An example of this index for various calculations at 6=1 

and 9 = 20°, 80° follows., In these calculations by the discrete- 

t 

angle method the nearest © i differed from © v by at least five degrees. 
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Another index used to test the calculations is reported 
in Section 5o22 0 

5.12 Form of the Equations for Calculation,, 

The expansion equations used in the method of integra- 
tion are as follows; 



m 


n/2 


sin k© 


7t/2 


sin k© 




/ f 2m (6,9) 


d© 


- ' ^ca< 5 ' 95 


d© 


m=l 


O 


cos j© 


O 


cos j© 






it/2 


sin k© 





- / y (6,n/2)sin© d© 

o cos j© 

k=l ^2f3t 0# *fk 

j=0, 1 ,2, o • . t (j— 1) 

m=k + j . 
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for a -2m ; 
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*/ 2 


sin k© 




m/2 




sin 


k© 


z W 6) 


/ f- (6,©) 

J 2m * 




d© 


= 1 Y sa 


(6,©) 




d© 


m=l 


o 


cos j© 




o 




cos 




tc/2 






sin k© 








- / 


V ( 

9 sa' 


6, 


n/2) sin© 




d9 






0 


OC*. 






cos j© 







k— 1 , 2 1 3 « • • • *k 
0 »l»2,«..,(j — 1 ) 

m=k + j » 



All integrations in these equations were reduced to the form 
of rapidly converging series or to algebraic results, before doing 
machine calculations. The forms used are given in Appendix A-3; the 
integrals on the right-hand side are not widely published. 



The expansion equations used in the discrete-angle method are 
given in Section 5oll, above. 



5° 2 The Principal Dependent Variables. 



5.21 The Results for A(6), B(6). 



The dependent variables A(6) and B(6) are calculated from: 



A(6) = ^ p 2m (6)Y 2ma (6 » u/2) + Y ca (6 » K/2) 

m-1 



m / , Nim-1 

“ 6 z ♦ ,lsin6 

m-1 



m 

B(6). = I q 2ffl (&) Y 2ma (6,7c/2) + Y sa (6,ll/2) 
m=l 



=6 f q_ (6) - [§• + Si(6)] cos (6) + Ci(&) sin (6) . 

jT-t 2m-l 2m 2 

m=i. 



The substitution for y (b 9 n/2) employs the identity 
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CO , ~x 

/ — Ij — — ~2 ^ “ Si (6)] cos { 6) + Ci(6) sin ( 6 ) , 

o 6 + x 

where Si(6) and Ci(6) are the sine and cosine integrals and are 
available in tabulated or series form, 

2 2 

The quantities A(6) and B(b) in the form A + B are use- 
ful for interpolation of the added~mass and damping coefficients, as 
pointed cut by Ursell [19571 , The phase angle of the motion 
e = tan~' B(b)/A(6) is smoothly varying and nearly linear in 6. The 
results of this calculation and that of Ursell [19493 for A(6), B(6) 
are given graphically in Figures 1 and 2, 

5,22 The Quantities M(6,0), N(6,©), and M q (6), N q (6)o 

The dependent variables M(6,9) and N(6,©) are required for 
the calculation of the hydrodynamic pressure distribution. They may 
be calculated from their definitions in Section 3° 4, which in this 
case become; 

= Z. q 2B (6) <P 2mi (6-S) ♦ <P sa <6-e> 

m=l 

0 

m 6 

~ 7 [cos 2si©+ cos (2m-l)G] + 

23L 

-6 cos O^Q + j sin n @) s ± n (6 sin q) + (^ + log 6 + 

n~x 

n 

+ 7 eoe n0) cos ( 6 sin 0) ] 9 
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FIGURE I. THE FUNCTIONS A(8) and B(8) 



tarr'(B(8)/A(8) 



46 




8 



FIGURE 2. THE FUNCTION tarr'(B(8)/A(8)) 
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N(5,©) = I P 2m (6) <P 2 ma (6 ’ 9) + ^ca 
m=l 

m 

= Z. P 2 m^ cos + gm l cos (2m-l)©] + n£~^ COS ®cos (6 sin ©) . 
m=l 

Here we have substituted the definitions of <p (6,©) and (p (6,©) 

CO. SB 

from Section 4*2 and a form suitable for numerical calculation of 

cp (6 f 9) from Appendix A-3*4* The open-ended summations over the in- 
sa 

dex n in M(6,9) were automatically run to "convergence" on the dig- 
ital computer, where "convergence" means the least significant digit in 
the floating-point calculation. 

The M(6,9) and N(6,©) were calculated for each 6 at inter- 
vals on the cylinder surface never larger than © = 0 O (10°)90°. 

The values of M(6,©) and N(6,0) available at © = 0 O (10°)90° 
were used to calculate M q (&) and N q ( 6) by a numerical quadrature. 

The defining equations are, from Section 3°8, 

n/2 

M (6) = ~ / M( 6,9) W(©)d© 

O i U 

o 

rc/2 

N (6) » ~ / N( 6 , 9) W(©)d© . 
o u 

o 

In this case G=1 and W(9) = cos 9. It is possible in the case of. 

the circular cylinder in infinitely deep fluid, as shown by Ursell [1949], 

to substitute the defining equations for M(6,9), N(6,9) and to simplify 

before calculating M , N . Here, however, we are concerned with the 

o o 

values of M(6,9), N(6,9) which other calculators bypassed. Using the 
calculated values of M(6,9), N(6,9) in a numerical calculation of 
M q (6), N q ( 6 ) gives an important check since (Section 3. 92) 
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M A - N B * n 2 /2 o 
o o 

The procedure employed was to perform a Lagrangian interpolation 
based on the calculated values of M(6,Q), N(6,©) at © = 0°(10°)90° 
and then to use Simpson’s Rule for approximate integration of 
M(6,©) W(©) and N(6,6) W(©) with a convenient interval of ©„ 

The calculations were tested at each 6 by computing 

M (6)A(6) - N (6) B(6) and observing the percentage deviation of 

o o t 

2 

the calculated value from n '2 Q Curves plotted through error index 
points are shown in Figure 3 as a function of the frequency parameter 
for sample calculations with six and nine expansion coefficients (each) 
p^Cd)* ^2n/^° ^t higher values of 6, seventeen expansion coeffi- 
cients were calculated by the discrete~angle methods 

5o3 Hydrodynamic Pressure Distribution on the Cylinder Surface 0 

The hydrodynamic pressure at the location 9 on the cylinder 
surface in units of the hydrostatic fluctuation is (Section 3o7) 

ELSuSI , 6 m;6,6;bW ♦ N'6,9:.A(6) (unit acceleration) + 

/° gh 0 A(6) 2 ♦ B(6) 2 

+ , M( 6,Q)A( 6) - N(6,|)J(61 (unit velocity) . 

AC 6) 2 + B(6) 2 

This is a simple, but repetitive, algebraic calculation using the 
principal dependent variables,, The pressure fluctuation in phase 
with acceleration, 

c MB 4- NA 



and the 



pressure fluctuation in phase with the velocity, 
MA - NB 



*v = 6 72 



A + B 
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S 



FIGURE 3. CALCULATION ERROR INDEX vs FREQUENCY. 

(M A - N B) - ir 2 /2 
Index n ^ 11 x 100% 

71/2 

( 1) Six expansion coefficients. Method of integration. 

(2) Six expansion coefficients. Discrete angle method. 

(3) Nine expansion coefficients. Discrete angle method. 
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were calculated for gives 6 and © - 0 C (10°)90 C ® 

Curves through calculated points of p a (6,©) and p v (6,©) are 
shown for example and Table 1 lists calculated values. 

These computed values are used for a prediction of the total 
pressure fluctuation on the cylinder used in the experimental work re- 
ported in Part C. 

It can be shown by deriving the velocity potential appropriate 
to the case of large 6 that the angular dependence of the hydrody- 
namic pressure approaches cos©, that is, for large 6 



Therefore 






6 cos © (unit acceleration) . 




MB + NA 



A 2 + 



B 



COS 0 



as 6 ■» oo . 



The sample calculations illustrate this since it is seen from Figure 6 
p^ tends to zero and from Figure 5 that p a /& tends to cos© for 
large values of 6® The same result holds for any cylinder of the el- 
liptic family® 



5ok The Added-Mass and Damping Coefficients® 

The ratio of the added mass to the mass of fluid displaced by 
the cylinder (Section 3<>9l) 

, M B + N A 
m 4 o o 

2 " !t .2 2 

ma yO/ 2 A -> B 

and the damping-force coefficient (Section 3® 92) 

F v _ 2n 

2 “ .2 „2 

na p/2 A + B 
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FIGURE 4. p (6,0) FOR THE CIRCULAR CYLINDER, 
a 
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FIGURE 5. THE TREND OF p^/8 FOR LARGE 6. 



53 




FIGURES, p (6,0) FOR THE CIRCULAR CYLINDER, 
v 
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TABLE I- A p (S.G) FOR THE CIRCULAR CYLINDER, 
a 
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were previously reported by Ursell [1949] for certain 0 <_ b < 3ti/2o 

His values are confirmed and the range of frequency parameter is ex- 

tendedo The values of M (6), N (6) used here were calculated from 

o o 

the values of M(6 ,©), N(6 ,©) used in the hydrodynamic pressure cal- 
culation as discussed in Section 5°22 0 The results are shown graph- 
ically in Figure 7o 



The results for the damping coefficient can be converted to 
amplitude ratio by 



h A 2 . 

, w, _ 6 it , 2v •. 

'■h ' ~ 2 ' 2 d 2 ; 

o A + B 



Ursell [1955 9 1954] derived the asymptotic result 



n i. 
w 

h ~ 6 
o 



as 6 ->• qd 



He pointed out in [1957] that although this asymptotic form is correct, 
it may be hazardous to interpolate with this result since earlier theo- 
ries predicted distinct minimao He suggested extended calculations 
could help settle the matter» The present calculations for the ampli- 
tude ratio tend distinctly to the asymptotic value as given by Ursell 
and shown in Figure 8» 



COEFFICIENT 
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FIGURE 7. THE ADDED-MASS AND DAMPING 
COEFFICIENTS FOR THE CIRCULAR 
CYLINDER. 
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FIGURE 8. THE AMPLITUDE RATIO A FOR THE CIRCULAR 
CYLINDER. 
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Vic Elliptic and More General Cylinders in 
Infinitely Deep Fluid 

6.1 Description of the Method. 

Sample calculations were made for the hydrodynamic pressure 
distribution, added-mass and damping coefficients of elliptic and 
more general cylinders in infinitely deep fluid. An important ob- 
jective of the sample calculations, in addition to that of obtain- 
ing the hydrodynamic quantities of interest, was to test the flex- 
ibility of the program written for control of the automatic digital 
computer. The program was successfully tested over a wide range of 
cylindrical shapes generally encompassing those of naval architec- 
tural interest. Some of these sample calculations are reported 
here. It is intended that a more comprehensive atlas will be sepa- 
rately reported later. 

The calculations were done on the IBM 704 digital computer at 
the Computer Center, University of California, Berkeley. 

6.11 Numerical Procedures. 

The expansion coefficients p^CS), ^2m^^ were calcu- 
lated by the discrete-angle method described in Section 5.11. The 
expansion equations are given in Section 4.3. The forms of 
Y (6,9) and y (6,m/2) suitable to numerical calculation are 
given in Appendix A«=3<>4o The number of expansion equations m and 
the discrete angles were varied® The least number m in the 

calculations reported here was nine so that the hydrodynamic pres- 
sure could be calculated at ten degree intervals with expansion 
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coefficients known to satisfy the truncated expansion equations 
exactly at ten degree intervals,, More often in these calcula- 
tions m = 17 ° In. every case the expansion coefficients were 
tested at arbitrary 9 = ©^ ^ ©^ to provide the accuracy index 
described in Section 5°11<> 



The principal dependent variables A(6), B(6), M(6,©), and 
N(6 t 9) were calculated in the same manner as for the case of the 
circular cylinder but with appropriate definitions of Y2 ma (6»Ti/2) 
and ^2ma^ & ’®^ aS S iven in Section 4,1, The dependent variables 
M q ( 6 ), N q ( 6) were calculated by numerical quadrature from the cal- 
culated values of M(6,©), N(6,©) as described in Section 5«22 

2 

and the result was tested by comparing M q A - N q B with rt /2 a 



6,12 The Selection of Mapping Parameters a 2 n+ i° 

The number N + 1 of mapping parameters 
a 2n+l’ n= ® *1»2, , . o ,N (Sections 2,2, 3,2) has been left unspecified 
in the problem solution,. It is convenient here to classify the 
cylinders by the number N + 1, The naval architect is more ac- 
customed to classifying ship sections by other parameters so some 
of this nomenclature is introduced. 



The ratio of the half-breadth of the cylinder at the free 
surface to the depth of the cylinder at the centerline is designated 
H: 

H = x(n/2) / y(0) = b/d . 



The ratio of the cross-sectional area of the cylinder to the area 
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2x( Ti/2)y(0) = 2bd, 



A 

is designated S and called the "area coefficient." 

The one-parameter cylinders generated by assigning N = 0 
suitable values to a^ 0 are members of the elliptic family. 

The special case a^ = 0 is the circular cylinder. The ratio 
H = (1 + a^) / (1 - a^) is uniquely determined. The area coeffi- 
cient for any elliptic cylinder is constant; S = n/4. Thus in 
naval architectural terms the one-parameter family of cylinders 
offers a choice of beam-draft ratio but constant area-coefficients. 



(2n+l)2 



/ <£ du, 



O _ 2L L_ ~ ‘ 
2 G 2 Po 



N 

L 

n=0 



The two-parameter cylinders generated by assigning N = 1 
and suitable values to a^ ^ 0, a^ ^ 0, to provide more-or-less 
ship-like cross-sections are commonly called Lewis forms, the ratio 
H and the area coefficient S, 



H = 



1 + + a^ 

1 — a^ 



S 



n 2 2 

71 1 ~ a ! 3a 3 
Tl ? 2 ’ 

4 (1 + a 3 r - a^ 



are determined by a^, a^ and conversely for ship-like forms. Thus 
the two-parameter family of Lewis forms permits some selection of 
area coefficient at fixed beam-draft ratio. 



Still greater flexibility in selecting a geometric form with 



given beam-draft ratio and area coefficient is possible in the 
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three-parameter family. Here N = 2 and appropriate values are 
assigned a^, a^, a^. L. Landweber and M. Macagno [19593 intro- 
duced a systematic treatment of these forms and present a conven- 
ient method for finding the asymptotic value of the added mass 
m^ (Section k,k). The introduction of the third parameter specif- 
ically permits a variation in the vertical distribution of the cross- 
section area at given beam-draft ratio and area-coefficient. This, 
in naval-architectural terms, is a two-dimensional interpretation of 
the vertical prismatic coefficient. Ship-like forms with higher area 
coefficient are also possible. 

Prohaska [19^73 considered two special cases from the four 
parameter family, N = 3* He selected non-zero values for a^, a^ 
and for a^, a^. The examples be illustrates of the case a^, a^ 
not zero are cusped, while several in the group a^, a ^ not zero 
may be called ship -like. 

The program for the control of the automatic computations 
permits the selection of any foreseeable number of coefficients. 
Coefficients are stated consecutively up to the highest non-zero 
value and the summations over N terminate at that number. Many 
sample calculations have been done with various combinations of 
coefficients through a^. The results of sample calculations of 
examples from the one-, two-, and three-parameter families are pre- 
sented here. 

6.2 Results for Examples from the One-Parameter (Elliptic) 

Family of Cylinders. 



The sample calculations for the circular cylinder given in 
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Chapter 5 are a special case of the one-parameter family of ellip- 
tic cylinders. The results of calculations for elliptic cylinders 
with H = 1/5 and H = 3/2 are presented in the following figures. 

The hydrodynamic pressure fluctuation in phase with the ac- 
celeration, p , and that in phase with the velocity, p (Sec- 
tion 5»3) are shown for typical values of © measured in the ref- 
erence plane as a function 6 = Kx(rc/2) = Kb. 

For large 6 m , 

— = t!L^g)^|.). _ cose 

6 A( 6; + B(6) 2 

for all cylinders of this family (Section 5«3). This fact is illus- 
trated by observing the slope of the p (6.0) lines and is particu- 
larly evident in Figure 10. 

As pointed out in Section 4.4 the free-surface correction 
k^ = m/m^ is numerically equal to the ratio of the added mass m to 
the mass of the fluid displaced by a circular cylinder of the same 
beam in the case of elliptic cylinders. Therefore curves of k^(6) 
show the frequency-dependence of the inertia-coefficient of Lewis 
[1929] or the added-mass coefficient of Landweber and Macagno [19573. 

Curves of the ratio of waveheight far from the cylinder to the 
amplitude of vertical oscillation, 

A = h /h 
w o 

are shown in Figure 14. The damping coefficient may be calculated 
from the amplitude ratio (Section 5.4), 
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FIGURE 9. 



p (6,0) FOR THE ELLIPSE H=3/2. 
a 
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FIGURE 10. p (6,9) FOR THE ELLIPSE H = 1/5. 
a 
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FIGURE ll» p (8,0) FOR THE ELLIPSE H= 3/2. 

v 




8 

FIGURE 12. p ( 8,0) FOR THE ELLIPSE H= 1/5. 

V 
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FIGURE 13. K FOR THE ELLIPTIC CYLINDERS H s l/5, 3/2. 
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FIGURE 14. A FOR THREE ELLIPTIC CYLINDERS. 
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Tasai [19593 has previously reported calculations for the free- 
surface correction factor k^ and the amplitude ratio A for the two 
ellipses reported here over a limited range of the frequency param- 
eter . His method of calculating six expansion coefficients satisfy- 
ing eight expansion equations in the least-squares sense was reported 
in Section 5.11. His procedure for evaluating the several integrals 
was not reported. The trends of all his calculations are confirmed 
and differences within his range of calculated values (determined 
graphically from his figures) are almost certainly due to differences 
in numerical procedures, plus the graphical dispersion. 

These curves give insight into the frequency dependence of the 
various hydrodynamic quantities for forms of various beam-draft ratio 
and constant area-coefficient. The vertical distribution of area as 
measured by the fractional depth of the section centroid or the frac- 
tional depth of the radius of gyration about the free-surface (x-axis) 
is also constant. 

6.3 Results for Examples From the Two-Parameter (Lewis Form) 

Family of Cylinders. 

The Lewis forms have, for given beam-draft ratio, a range of 
values of area-coefficient. The vertical distribution of area of the 
forms of fixed beam-draft ratio also changes ( cf . the three-parameter 
family of the next section). This one added mapping parameter fortui- 
tously provides sections that are more full than the ellipse and gen- 
erally accepted as ship-like. The upper limit of fullness obtainable 
is somewhat arbitrary. Lewis [19293 preferred to judge the forms by 
the practiced eye of a skilled naval architect; Landweber and 
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Macagno [1957^ formalized this decision. The lower limit of fullness 
is probably also arbitrary since the forms tend to be not aestheti- 
cally ship-like and if fine enough, are cusped at (3 = 0 for H <_ 1 
and at p = n/2 for H > 1, 

The sample calculations made to date suggest that the behavior 
of the various hydrodynamic quantities does not change abruptly as 
the values of a^, a^ encompass the reasonably expected range. 

The results of two of the calculations for this family are 
shown graphically. Both forms have H = 1, the same as the circular 
cylinder. One form is full and ship-like. The other is the section 
of minimum area and is cusped. 

Figures 15 and 16 show curves through calculated points of 

p (6,0). The asymptotic behavior of p /6 for large 6 for this 
a a 

family is 



p (1+a ) cos© - a cos 3© 

— - ~ as 6 oo . 

6 C l+a n + a^ ) 

1 5 

This behavior, more complex than simply cos 9 as in the one-parameter 
family, is shown by the slopes of the curves in the two figures* 

Figure 19 shows curves of the free-surface correction factor, 
k; — m/m . It is recalled that for the more general cylinders the 
free-surface correction factor k^ is no longer numerically equal 
to the added-mass or inertia-coefficient* The calculation of Tasai 
C1959] for k^ of the full-form is shown over his range 
0*4 < 6 < 2.6. 

The damping coefficients for the full form and the cusped form 



are shown in Figure 20 
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FIGURE 15. p (8,0) FOR THE FULL SECTION, H= I. 
a 
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FIGURE 16. p (6,0) FOR THE CUSPED SECTION, H= I . 
a 



6 
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FIGURE 17. p v (6,0) FOR THE FULL SECTION, H=l. 




6 

FIGURE 18. p (6,0) FOR THE CUSPED SECTION, H=l. 
v 
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FIGURE 19. k FOR THREE CYLINDERS, H»l. 
4 



DAMPING COEFFICIENT 
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FIGURE 20. 



THE DAMPING COEFFICIENT ' " 

A 2 + B2 

FOR THREE CYLINDERS, H = l. 



74 



6.4 Results for Examples From a Three-Parameter (Landweber-Macagno) 

Family of Cylinders. 

The use of three parameters a^ a^, a in the mapping func- 
tion permits some independent variation of the vertical distribution 
of area for the ship-like section of given beam-draft ratio and area 
coefficient. L. Landweber and M. Macagno [19593 chose to measure the 
vertical distribution of area in terms of the parameter q defined 
as the moment of inertia of the actual section about x-axis made non- 

z z 

dimensional by the product x(u/2) y(0) = bd . Another useful char- 

acteristic of the three-parameter family is that it provides ship- 
like sections of larger area-coefficient than is possible in the two- 
parameter family. 

Examples are shown which illustrate the behavior of the hydro- 
dynamic quantities for a bulb-like section and a section of maximum 
fullness. The sections are shown in Figure 21. The full section is 
an example of sections possible in the three-parameter family that 
are ship-like and have area coefficients larger than is possible in 
Lewis forms. The area coefficient for the full section is S = 1.012. 
This section was selected from Figure 4 of the Landweber and Macagno 
reference cited above. The bulb-like section has an area-coefficient 
of about 0.6 and it is a member of the special group a^ =0; a^, a,- 
not zero from which Prohaska [194?] selected examples. 

The amplitude of the pressure fluctuation in phase with the ac- 
celeration at locations near the vertical centerline of the full sec- 
tion increases much more rapidly with increased 6 = Kb than for the 
more rounded elliptic or bulbous forms. 
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FIGURE 21. 


THREE CYLINDERS, 


H s 1/5 • 
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FIGURE 22. p Q (6,e). FIGURE 23. p y (6,0). 

p (6,0) AND p (6,0) FOR THE FULL FORM, H= 1/5 . 
a v 
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The asymptotic form of p /6 for large 6 for the bulbous 

cX 

form is shown in Figure 2k. 

In Figure 27i showing the free surface correction factor 
and the amplitude ratio A, data on the elliptic form with the same 
beam-draft ratio, previously given in Section 6*2, are reproduced 
for reference. The amplitude ratio for the bulbous shape indicates 
a relatively high damping coefficient* the full section (larger 
vertical prismatic coefficient) is relatively less damped in heaving 
oscillation. 
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FIGURE 24. p (6,©) FOR THE BULBOUS SECTION, H=l/5. 
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FIGURE 25. THE TREND OF p /6 FOR LARGE 6, 

a 

BULBOUS SECTION, H = 1/5. 
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FIGURE 26. p (6,©) FOR THE BULBOUS SECTION, H= 1/5. 
v 
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FIGURE 27. k 4 AND A FOR THREE CYLINDERS, H* 1/5. 
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VI I . Experiment and Results 
7.1 Purpose and Direction,, 

F. Urseil, who obtained the first rigorous solution for the 
circular cylinder in fluid ©f infinite depth [1949] also said clearly 
"experiments are needed" and has reiterated this even in recent times: 

"It would be highly desirable to compare these calculations 

with suitable experiments,," [19573 . 

The response has not been overwhelming,, The earlier experiments 
reported in Chapter I are not strictly comparable, and in later 
years, three-dimensional experiments (Golovato [1957], Gerritsma [1957]) 
have received more attention. Apparently, the only modern experiment 
simulating two-dimensional behavior is that of Yu [i960], who mea- 
sured waveheight ratio and found satisfactory agreement while the pres- 
ent experimental apparatus was under construction, Yu reported instru- 
mentation difficulty in measuring the vertical force acting on the 
cylinder. 

Pressure fluctuation measurement was selected as the primary 
objective of the experiment to be discussed. The total pressure 
fluctuation signalled by a suitable transducer is a distinctive func- 
tion of the location 0 on the cylinder, as well as the frequency 6, 
according to the preceding theory,, It is considered unlikely that 
reasonable agreement could be expected between experimental results 
and the distinctive prediction unless the analytic model accurately 
represents the essentials cf physical reality. Vertical force acting 
on the cylinder was selected as a secondary objective because, first, 
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this is an integrated result of a predictable primary physical quan- 
tity, Second, it appears not feasible to avoid encumbering the mea- 
surement of force by ordinary body inertia reaction, which is far 
from negligible, yet not of concern to this theory. 

An additional reason for placing emphasis on the pressure 
distribution measurement is the expectation that such investigation 
will help to provide insight into similar problems of naval archi- 
tectural interest. 

The following measurements were obtained with a cylinder of 
10-inch radius in forced vertical harmonic motion of small ampli- 
tude. Details of the experimental arrangement and simulation of two- 
dimensional conditions are recorded in later sections. Over a cer- 
tain range of oscillation frequency, 0,2 < f < 0,6 cycles per second, 
the particular arrangement used simulated two-dimensional, conditions 
less successfully and this caused consistent small deviations from 
theory in this range. The experiments were conducted at the Ship 
Model Towing Tank, Richmond Field Station, University of California, 

7.2 Results of Measurements of the Pressure Fluctuation, 

The total pressure fluctuation at a point on the surface of 
a horizontal circular cylinder given forced vertical oscillations in 
an ideal fluid of infinite depth is predicted with the aid of the 
calculations reported in Section 5 ° 3 » The total pressure fluctuation, 
in units of the hydrostatic fluctuation for a forced motion y = h^ sin cot, 
is 

- s [l-p (6,©)] sin cot + p„(6,©) coscot . 
p gh, ' 
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Here p (6,9) and p..(6,9) are the acceleratioii- and velocity- 

cl v 

phase components of the hydrodynamic pressure given for this case 
in Section 5o'j>» The amplitude of the total pressure fluctuation is 



p( 6 , ©) 



Eh& s £l 

pg h 0 



[(!• 



“ p a } 




o 



The predicted results are shown by solid lines through cal- 
culated points in Figure 28 for 9 = 0°, 40°, and 60'"' o Experi- 
mental points for measurements at the same locations are shown on 
the same figure 0 These points are (amplitude of the recorded 
pressure fluctuation)/( amplitude of the recorded pressure fluctua- 
tion for hydrostatic displacement) „ The distinctive trends, for 
each ©, of the predicted and measured results are generally con- 
fluent o An observed departure from two-dimensional conditions in 
the experiment is the pz*incipal cause for consistent small devia- 
tions from predicted results in the range 0,2 < f < 0,6 cycles/second. 
The trend toward higher-than-predicted results at © = 60 at the 
higher frequencies observed may be caused by viscous effects in the 
locations of higher tangential velocity. 



7o3 Results of Measurements of the Vertical Force 
Acting on the Cylinder, 

The total vertical force per unit length required to sustain 
steady-state motion y - sin cut of the test cylinder is the sum of 
the hydrodynamic force (Section 3o8), 

_ M B + N A ti 2,„ 

_ 2 o o " „ 2 m /2 • 

2 y + ojy , 

A + a A + B 
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FIGURE 28. THE CALCULATED AND MEASURED AMPLITUDE OF TOTAL PRESSURE 
FLUCTUATION. 
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and the hydrostatic force 2aog y and the inertia reaction of the 
mass of the cylinder. The mass of the cylinder in this experiment 
was made equal to that of the displaced fluid for neutral buoyancy 
at rest. The total force in units of the change in hydrostatic 
force is 



F(t) 

2pgah^ 



Cl - (l+k)6it/4] cos 0 )t ~ ~2 — sin cut 

A + 



k - 



4 

T£ 



M B + N A 
o o 




o 



The amplitude of the total force fluctuation is the square-root of 



[1 - (l+k)6n:/4] 2 + (6 tsA)3 2 . 

A" + B 

The added mass coefficient k and the damping force coefficient 
2 2 

2u/(A + B ) for this case are reported in Section 

The predicted result is shown by the solid line through cal- 
culated points in Figure 29- Points for measured values are shown. 
These points are (amplitude of the recorded total force signal) 
divided by (amplitude of the recorded signal for hydrostatic dis- 
placement). The general agreement is satisfactory. Again, the 
observed departure from two-dimensional conditions is the principal 
cause for small differences from predicted results in the range 
0.2 < f < 0.6 cycles/second. 

Measured results are generally higher than predicted in the 
neighborhood of the force minimum around f -- 0.9 epa. In this 
region 



NORMALIZED VERTICAL FORCE 
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FIGURE 29. CALCULATED AND MEASURED TOTAL 
VERTICAL FORCE. 
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1 - (l+k)6n/4 = 0 

and the total force is nearly that in phase with the velocity. Any 
additional energy dissipation (due to viscous effects, end-effects, 
etc.) would contribute additional force increments in phase with the 
velocity and not predicted by the theory. The measured values there- 
fore may well be expected to be higher than predicted and more dis- 
tinctly so when the much larger component in time-quadrature vanishes. 

7.4 The Cylinder and The Experiment Arrangement. 

7.4l The Cylinder and The Motion Generator. 

The cylinder was fabricated from a one-quarter inch thick 

plate of aluminum rolled to a circular submerged cross-section with 
10-inch outside radius and 12-inch vertical sides above water. The 
cylinder is 4l inches long. The ends are closed by end-plates of 
three-quarter inch aluminum plate. Two heavy steel athwartship 
frames, machined to semicircular shape, together with the end frames, 
assure semi-circular shape, contribute rigidity, and anchor parts of 
two force dynamometers. The model plus the parts of the dynamometer 
in the model was weighted to neutral buoyancy at the desired water- 
line. The high vertical walls anchor dynamometer parts, contribute 
torsional rigidity and may be used in other experiments with deeper 
submergence. These high weights cause static instability but the 
overturning moment is small and of little consequence since the model 
is secured by means of vertical -motion dynamometers to the rigid 
motion generator. 

The vertical-motion force dynamometers consist of an inter- 
changeable set of beam "springs" each 7-1/2 inches long, 2-inches 
wide and 5/32 or 3/l6-inches thick for choice of spring-constant. 
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FIGURE 30. 



THE 



MODEL 
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FIGURE 31. THE VERTICAL MOTION DYNAMOMETER 
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Ah array ®f six ®ucfc springs connects each ®£ two mottost-generator re- 
ciprocating rods to it© steel-frame aad end-plate in th® cylinder la 
as array strongly preferential to vertical relative motion* The rela- 
tive displacement between the model and the reciprocating rod, a mea— 
sure of the vertical force acting, was sensed by a Statham Instruments 
Corporation Model G— IQ— 3 unbonded strain-gauge displacement transducer* 
Relative vertical motions of about 0*002 inch were permitted and this 
is considered a negligible uncertainty in th® stroke of the motion 
ganeyator 0 The ©aaplituda of the stroke was varied from ©ae-half to 
on# inch and had n© observable effect on the recorded. Measurements* 
Almost all of the pressure measurements were made with one-half inch 
stroke* 

The motion generator is an elaborate machine designed and 
build for th® Ship Model Towing Tank* A 3/4=horsepow«r constant^ 
speed motor,, through a continuously -variable speed r®du<f®r, drives 
a pair of mechanisms patterned after th© "’scotch-yoke 80 which trans- 
late eoastant-speed rotation to uniform simple-harmonic reciproca- 
tion of two stool rods* The scotch-yoke, at th® expense of elabora- 
tion , has th® advantage of providing almost exact sinusoidal motion* 

The motion was continuously monitored by raacha&ically-driv®® poten- 
tiometers or the mechanism providing electrical signals which,, be- 
cause of orientation, were sine and cosine signals* 

The pressure gauges were installed ia the shell of the cylin- 
der at various angular locations at mid-length* Th® gauges were un- 
bonded strain-gauge types with sensitive diaphrams about 3/4-inch 
diameter* The diaphragms therefore subtend (Oo75)/XO m 0*075 radian, 
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about 4.3 degrees and it is understood results measured "at 40 degrees" 
admit the diaphragm extended about +_ 2 degrees* The Statham Instru- 
ments Corporation Type PM-233“TC pressure transducers have calibra- 
tion factors close to 23 millivolts/psi with strain-gauge bridge exci- 
tation of 5 volts* The full-scale range is +_ 0.5 psi differential. 

The differential feature was used to provide rear-side pressurization 
(by air and static manometer) to balance the hydrostatic pressure pre- 
load on the sensitive diaphragm* This contributes to maximum initial 
sensitivity and linearity? the response of diaphragm-type transducers 
was studied extensively and is reported in the next chapter. The dy- 
namic response of these gauges is uniform from hydrostatic up through 
frequencies well above those reported here 0 

7.42 The Simulation of Two-Dimensional Conditions. 

The experiments were conducted in the Ship Model Towing Tank, 
Richmond Field Station, University of California* The tank is 200 
feet long, 8 feet wide and 6 feet deep* The cylinder and motion 
generator were installed mid-length and transversely between false 
walls installed in the tank to narrow the width at the model to 42 
inches. The model is 4l inches long so there was 1/2 inch clearance 
at the ends. At the frequencies observed, surface action and dis- 
turbances at the ends was not violent but was increasing with fre- 
quency. At frequencies above those reported here the surface wave 
action was not two-dimensional. 

The false walls at the model were 8 feet long, extended above 
the surface, and were solid to the bottom. The walls were 5-pound 
steel plate rigidly stiffened on the outside by welded flatbar 
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stiffeners , The para.i e f a wails at tne mod°l were joined to the 
regular tank walls by further a- : + ions at each end about- 8 feet long. 
Thus the model was in the 42 In oh spa's between false walls which were 
parallel for 8 feet and flared out on each side to 96 inch spacing. 

It is agreed that for the experiment reported here alone, a model 8 
feet long without special walls would be preferred*, but other work 
was also in mind. Nevertheless, with the exception discussed below, 
it appeared that two-dimensional conditions generally prevailed and 
regular progressive surface waves were generated. Since the motion- 
generator was over-powered and positively driven, steady-state condi- 
tions were approached in the neighborhood of the cylinder very rapidly. 

Under certain conditions, described below, a distinct departure 
from two-dimensional conditions was observed. The influence of this 
effect as reflected for example in the total force measurement was 
strongly dependent on the transition conditions at the ends of the 
parallel walls. Three transition conditions are shown in Figure 32, 

Arrangement (a) is the least abrupt transition from k2 inch 
spacing at the model to the 96 inch tank width. This is the arrange- 
ment used for the data recorded in the preceding sections. The con- 
sistent lower-than- theoretical experimental values in the range about 
0,2 < f < 0,6 cyclee/eecond are believed to be caused principally by 
this transition. The effect is much mere intense with more abrupt 
transitions. Two other arrangements were tried and are illustrated 
in the figure together with sketches shewing the trend of total, force 
measurements in the criti.-a.’ fr^cuen y range. The consistent pattern 
leads to the belief that parallel walls with a fitted model would 



94 






This side 
same as (a) 



ARRANGEMENT (b) 





ARRANGEMENT (c) 




FIGURE TOP VIEW OF ARRANGEMENTS, (b) and 

(c) accentuate departure from two dimensional 
conditions. 
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eliminate this cause of deviation of measured from predicted values. 
In all arrangements the surface wave-action, even in this critical 
range, appeared to be normal and two-dimensional in the near neigh- 
borhood of the model except for minor disturbances in the corners at 
the ends. 

7.4} Electrical Instrumentation. 

The unbonded strain-gauge displacement transducers in the 
force dynamometer and pressure transducers in the shell of the model 
(Section 7»4l) were connected to Brush Instruments Division, Clevite 
Corporation, Model RD-5612-00 carrier and recording amplifiers. The 
carrier frequency is nominally 2 kcps, the bridge excitation voltage 
normally available is less than maximum usable by the transducers 
but signal strength was sufficient and reliable. 

The electrical signals monitoring the motion (Section 7.41) 
were connected to the direct-current recording amplifier sections 
of other Brush RD-5612-00 units. 

All electrical signals were recorded on an 8-channel Brush 
Instruments Model RD-2684-50 linear chart recorder. Data were re- 
duced from the records manually; recorded amplitudes measured in the 
range 0.6 to 4 centimeters. 

The frequency of oscillation was computed from the recorded 
motion sinusoids and known chart-speeds and also compared with rev- 
olutions counted by a Berkeley Instruments Corporation "events-per- 
unit-time n counter. This instrument was given electrical pulses 
from cam-operated microswitches and the number of counts in a known 
time interval confirmed the oscillation frequency. 
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VIII. General Study of Instrumentation for Measurement of 
Pressure Fluctuations. 



8.1 General Characteristics of the Pressure Measurement. 



This chapter presents an analysis of pressure transducers 
and their use in measuring the hydrostatic and hydrodynamic pres- 
sure distribution on a body oscillating in a free surface. The 
material is part of a general study of instrumentation pertinent 
to the experimental investigation reported in the preceding chap- 
ter and is of interest in similar investigations. 



The predicted values for the experimental case may be used 
to fix ideas on the pressure fluctuation magnitudes and frequency 
range. The fluctuation is periodic at the model-motion frequency. 
Its amplitude is 



p(6,9) = 



p(6,9) 



/°eh 



= C(i-P a ) 2 + P v 2 ] 



1/2 



where p (6,9), p (6,9) are the acceleration- and velocity-phase 

9. V 

components calculated for the example in Section 5.3» Curves 
showing this result are given for various 9 in Figure 28. 



At low frequencies, 6 < 0.2, p(6,9) = 1 independent of 

location 9 . That is, the pressure fluctuation is approximately 
the hydrostatic change due to slow heaving. The surface is barely 
distorted, and the pressure change is nearly independent of Toca- 
tion. This is rigorously correct as 6 0. Since the displacement 

is assumed to be small with respect to body dimensions, this refer- 
ence pressure change is small. 
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The amplitude of the pressure fluctuation goes through a 
minimum at a characteristic frequency depending on the location 9. 

For example, the minimum is about 20% of the hydrostatic change at 
locations in the bottom sixty-degree sector and occurs at 1,5 < 6 < 5. 
Twenty percent of the expected hydrostatic pressure fluctuation is 
a useful index for the minimum expected pressure fluctuation. For 
the 20-inch-diameter cylinder used in the experimental part of this 
work, 20 percent of the one-half inch motion in water corresponds 
to about 0.004 psi. 

At higher frequencies, the magnitude of the pressure fluc- 
tuation increases. Ultimately, cavitation at the bottom center- 
line or viscous effects or short surface wavelength violate the 
assumptions preceding these predictions. 

The frequency range of interest in the experiments with a 
cylinder of 20 inch diameter was 0-2 cycles per second, that is, 
about 0 < 6 < 3» 

The low frequencies and preference to record hydrostatic 
pressure suggest a diaphragm-type transducer. Piezoelectric-type 
hydrophones are contra-indicated by this requirement because of 
their limited low frequency response. 

The low pressures and pressure fluctuations dictate high 

sensitivity. For ship-like cylinders of about 1-foot beam the 

minimum detectable pressure change should be about 0.001 psig (about 

2 

l/4o inch water or +110 DB ref 0.0002 dynes/cm ). This requirement 
is easily met by a piezoelectric transducer but the high sensitivity 
conflicts with the requirement for small diameter in the case of 
diaphragm types. 
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For experimental convenience, cylinder dimensions must be 
manageable and for maximum utility the transducers should be com- 
patible with the ship models ordinarily usedo 

8.2 Piezoelectric Transducers. 

8o21 Open Circuit Operation* the Piezoelectric Wafer. 

The piezoelectric transducer gives promise for this applica- 
tion because of its sensitivity o Its disadvantages are limited low- 
frequency response and sensitivity to local vibrations of its mount. 

To illustrate the operation of a piezoelectric transducer 
and its major disadvantage, consider the response of a disc to nor- 
mally incident pressure signals. The disc area is A, its thickness 
t. Assuming a material similar to barium titanate, the disc is po- 

i 

larized across the thickness direction, a. The rear face is blocked 
by a relatively large mass. 

This wafer of piezoelectric crystal, compressed in the direc- 
tion of its thickness t with a pressure p to a stress level 
p = F/A, generates a proportional electric field in the direction 
of polarization: 

E = kpai , 

E = V/t = kp . 

The potential difference across the wafer is 
V = (kt)p = chp 



c?" = pressure sensitivity, volts/microbar 
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The constant k in units of volts/cm thickness/microbar pressure 

-5 

is typically 1.4l x 10 for ceramic ,? B M Barium Titanate and 

-5 

2.4l x 10 for common piezoelectric materials called PZT-4 and 

2 

PZT-5 • (1 microbar = 1 dyne/cm ) 0 

The open circuit sensitivity of a piezoelectric transducer 
is therefore independent of diameter and is easily of the order of 
-110 db referred to 1 volt per microbar. Thus piezoelectric trans- 
ducers easily meet the required characteristics of high-sensitivity 
and small diameter. 

The required low frequency response requires closer examina- 
tion due to the loading effect on the crystal of the external cir- 
cuitry. To examine this, one uses the equivalent circuit of the 
crystal element. This consists of a voltage generator of zero in- 
ternal impedance in series with an impedance equal to the internal 
impedance of the crystal. The voltage output of the conceptual 
generator is equivalent to the crystal pressure sensitivity. The 
internal impedance of the crystal is considered constant and equal 
to that of the capacitor formed by the crystal of dielectric mate- 
rial between equal-potential conducting surfaces. 




e ( t ) = cTp(t) 



8.22 Low Frequency Response of a Piezoelectric Transducer 
and Amplifier. 

The input impedance of an amplifier connected to a piezo- 
electric transducer degrades the low-frequency response relative 
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to the open-circuit performance. This is seen by considering the 
amplifier input signal e^ : 




In the usual circuit-analysis technique we use 



s = joo 



E^(s) = transform of e^Ct) 
E(s) = transform of e(t) 

= cr P(s) 



E.(s) 



E( s ) 



s 

s + 1 /t 



tf P(s) 



s 

s + 1/t 



T = Rc o 



Consider a step-change in pressure: 

f . x fo t < 0 
P (t> =|a t > 0. 

The amplifier input signal is 

e i(t) = o' Ae~^ T (t >_ 0) o 

The open-circuit (unloaded) response of the crystal transducer 
would have been o' A • The penalty imposed by the external ampli- 
fier (load) is the exponential decay of the hydrostatic response 



with a time constant t 
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Response to a hydrostatic change in pressure. 



The steady-state response to a sinusoidal time-varying pres- 
sure decreases in amplitude and is increasing in phase error below 
the cutoff frequency. The cutoff frequency is 




Amplitude Response and Phase Error relative 
to Open-Circuit Response of a Piezoelectric 
Crystal for Sinusoidal Pressure Fluctuation 
of Frequency to . 
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Fractional Amplitude Response 

OF A PIEZOELECTRIC TRANSDUCER AND 
DC AMPLIFIER AS A FUNCTION OF 
FREQUENCY NEAR CUTOFF O) 0 = -5 



FIGURE. 37 . 
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These results show that the lew frequency response of the 



piezoelectric transducer plus amplifier is limited by the amplifier 
input impedance R and the crystal equivalent capacitance c, The 
amplifier input impedance is bounded by practical considerations of 
the associated thermal noise and by susceptibility to extraneous re- 
sponse due to the high impedance level© The crystal equivalent ca- 
pacitance c can be increased by increasing the wafer diameter, 
which conflicts with small size 5 or by decreasing wafer thickness 9 
which decreases sensitivity^ or by stacking thin wafers, which in- 
creases complexity of the transducer© 



ducers will be met with increasing difficulty as the desired low 
frequency response is extended* 

Some techniques are suggested to partly alleviate this limi- 
tation* 



8©23 Extending Low Frequency Response by Capacitive Loading. 
The low frequency response of piezoelectric transducer plus 



amplfier inay be extended at the expense of sensitivity. The ampli- 
fier input resistance R is shunted with a capacitance C© 



Therefore it can be expected the use of piezoelectric trans- 




c ; 




E ± (s) = E(s) 



c. __ 
C + c s 




s 




x = P.(C + c) 
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Comparing this result with that of the amplifier alone, 

(i) The apparent sensitivity is reduced by a factor 
c/(C+c) £ loO 

(ii) The hydrostatic time-decay factor is increased by 
a factor (C+c)/c© This is equivalent to reducing 
the low-frequency cutcff frequency by the factor 
c/(C+c) <_ 1.0 

With this technique the relative response is kept uniform 
to a new lower cutoff-frequency at the expense of being at a lower 
absolute level* The degree to which this can be pursued thus de- 
pends on the original crystal sensitivity and the magnitude of 
pressure fluctuations which are to produce measureable signals with- 
in the absolute level capability of the available amplifier© That 
is, the signal to noise ratio at the reduced apparent sensitivity 
will limit an attempt to indefinitely extend the low-frequency re- 
sponse by this technique© 

8©24 Extending Low-Frequency Response by Feedback 
Amplifier Technique 0 

Manipulations of the low-frequency cutoff and of the appar- 
ent sensitivity, interdependent in the preceding technique, can be 
largely separated at the expense of increased electronic complexity© 
An operational amplifier of standard design is introduced with ca- 
pacitive feedback© The amplifier will have a large gain, -A, with 
negative sense meaning phase reversal© The amplifier is realisti- 
cally assumed to have constant gain over frequencies of interest 
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and to have negligible output impedance# The input impedance cannot 
be neglected and as before is R 0 This amplifier is to be considered 
as a preamplifier or buffer for operational purposes between the pie- 
zoelectric transducer and the usual recording system© 




By the usual circuit-analysis methods and without further 
assumptions. 



E q (s) = - E(s) 




r 

\ 



1 + b 



s 



S Hr 



1 



T 



tf 



= - <s 






PCs) 



S + — 



T" = R[c + C (1+A)] . 



Comparing this result with those preceding it is seen the sensitivity 
and time constant are changed by functions of C and A. If the 
feedback element C is reduced to zero the relations are simply 
those for an amplifier following an amplifier-loaded piezoelectric 
transducer. If C is not zero and the realistic assumption is 



made that 
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| A | >> 1 

then the relations are simplified. Compared with the amplifier- 
loaded case, 

(i) the sensitivity is modified by the factor c/C 

(ii) the time constant is increased by the factor 
C(l+A)/c % CA/c . 

This is equivalent to reducing the low-frequency 
cutoff by a factor c/CA . 

It is striking to observe then that if C is made approximately 
equal to c, then 

(i) the sensitivity is equal to the open circuit 
sensitivity. 

(ii) the low-frequency cutoff is reduced by the factor 
1/A, where A is a large number. 

It appears this should be a useful technique for extending low- 
frequency response with nearly independent control of the appar- 
ent sensitivity and the extent of cutoff reduction. Practical 
limitation on the extent of this manipulation is not clear. 

8.3 The Pressure Sensitivity of Diaphragm-Type Transducers. 
8.31 Definition of Sensitivity. 

The transducers considered consist of two elements. 

(i) the mechanical pressure receiver. Primarily this 
will be a diaphragm but may be a tube, bellows or 
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C/c 

Fractional response, as a function 
of (added / internal) capacitance 



FIGURE 34 ‘ 



capsule which deflects under pressure and seals the 
device against fluid ingress. 

(ii) the deflection detection schema. This converts the 
pressure-receiver (diaphragm) deflection to a signal 
in useful form. In order to take advantage of elec- 
tronic (linear) amplification, electrical schemes 
are of most interest in devices of high sensitivity 
and resolution. Examples may be found using various 
physical principles: variable resistance, capacitance, 
reluctance, light flux 9 etc. 

Pressure variations on the pressure receivers of several different 
types of transducers change the internal energy in the detection 
device a negligible amount compared with the change in elastic strain 
energy in the pressure receiver. Limiting consideration to such 
stiffness-controlled transducers, it is possible to analyze the two 
elements of the transducer separately, and the transducer pressure 
sensitivity is the product of the sensitivities of each element: 



6 = 6(p), the diaphragm deflection at the center; 

e = e (6), the transducer signal output; 
o o 

de ,, de 
o _ do_ 

dp "dp d6 

The effectiveness of the pressure measurement system cannot be 
judged on transducer sensitivity de^/dp alone because this signal 
will be contaminated by a noise signal. Also, a response index will 
be required characterizing the least readable amplitude fluctuation 
of the indicating needle or chart line. Before attacking these 
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problems it is fundamental to study the deflection-pressure rela- 
tion 6 = 6 (p). This section will consider this relation and the 
pressure sensitivity, 



The theoretical maximum pressure sensitivity of a diaphragm 
is generally considered to be the initial sensitivity of a flat 
circular plate under uniform pressure load on one side producing 
infinitesimal inextensional strain,, 

By inextensional strain is meant that the plate neutral 
plane is not extended,. The internal strain energy is greater and 
the plate is stiffer if extensional strains exist* Inextensional 
conditions are consistent with Kirchhoff plate theory. 

The common and useful practice of diaphragm corrugations 
invariably introduces locations of extensional strain. Corrugated 
diaphragms realize less sensitivity than the theoretical initial 
sensitivity of the equal flat disc in Kirchhoff plate theory. 

The Kirchhoff theory, neglecting stretch and shear of the 
plate neutral plane, yields the following result for the deflec- 
tion 6 at the center of a clamped-edge flat disc of diameter D 
uniformly loaded by a pressure p on one side: 




8.32 Small Deflections of Clamped~Edge Flat Circular 
Plates 5 Kirchhoff Theory. 




where 







/ 
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2 3 ? 

F = a plate modulus E/(l-ji ), not Et /12(l-ji ) , 
t = thickness. 

The deflection is linear in pressure^ the pressure sensi- 
tivity is independent of the deflection. 

Non-dimensional expressions for deflection and for pres- 
sure are: 




These relations are illustrated graphically in Figure 35 
and 36. These representations are subject to the assumptions of 
their theory. Thus it is clearly not correct that 6/t increases 
linearly indefinitely and it is not possible to indefinitely in- 
crease sensitivity cTF/D by decreasing t/D o 

We will consider several cases illustrating that the pres- 
sure sensitivity actually achieved by a thin diaphragm will prob- 
ably be less than this theoretical maximum. 

8.33 Pressure Sensitivity of Flat Circular Plates with 
Deflection not Infinitesimal. 

For practical purposes it is not always possible to limit 
diaphragm deflections to deflections small in comparison with dia- 
phragm thicknesso Such factors as hydrostatic preload or practical 
deflection detection schemes lead to larger deflections and, as we 
will see, possibly decreased sensitivity. 



The exact solution of the equations of equilibrium of a flat 



di CO 



ill 




Z 4 6 <3 lo 

tsr 



Nondimensional Deflection of 

OF FLAT CIRCULAR DIAPHRAGMS 
WITH CLAMPED EDGES AND NO 
INITIAL TENSION 

(1) KIRCHHOFF 

(2) WAV 

(}) H E.NCKY 



FIGUPL 5^. 
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NON-DIMENSIONAL PRESSURE SENSITTVITY 
OF A FLAT CIRCULAR PLATE 



FIGURE 36 
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circular disc uniformly loaded on one side and with clamped edges 
and no initial tension is due to S. Way [1934]. This solution is 
available in tabular or graphic form but not in algebraic or series 
form. This solution is also illustrated in Figure 35. 

Approximate solutions are due to A. Nadai [ 1925 ], 

S. Timoshenko [1940], and A. Griffith [1927]. Another approximate 
solution is reported by E. Waters in the discussion accompanying 
the solution by Way [1934], These approximate solutions are ex- 
pressed in algebraic form. The deflection-pressure function is non- 
linear and additional definitions of sensitivity are useful. 

Pressure sensitivity is d&/dp as before and may be called 
the instantaneous pressure sensitivity. Initial sensitivity is de- 
fined as the limit of the pressure sensitivity for infinitesimal 
deflections: 



O'. 

1 



lim ^ 
p.o dp 



Average sensitivity at a given deflection and pressure is 



o' - 




The ideal flat-disc case has constant pressure sensitivity 

d6/dp equal to the initial sensitivity d6/dp] and average sen- 

p=o 

sitivity 6(p Q )/p o . 

The graphical presentation of the large deflection case sug- 
gests that the initial sensitivity is equal to that in the ideal case. 
The pressure sensitivity decreases for increasing deflection. The 
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average sensitivity at a given pressure, from algebraic expressions 
similar to those of the approximate solutions of Timoshenko, Nadai, 
Griffith or Waters is 




(|) 3 0.545 ( 



1 




[the constant 0<>545 may vary about +_ 10 percent depending on Poisson’s 



average and agrees with Griffith for y. = O. 29 ]. This is illustrated 
in Figure 36 for two assumed values of 



It is seen the average sensitivity may well be markedly lower than the 
theoretical maximum. The instantaneous pressure sensitivity is always 
less than the average sensitivity at finite deflections. 

Corrugated diaphragms such as those of the National Bureau of 
Standards series realize average sensitivities somewhere between the 
ideal maximum and the value of the average sensitivity for 6 (p q )D = 1/100. 

This limitation of sensitivity to values less than maximum 
theoretical applies when deflections are approaching 0.4t or more. 

High sensitivity requires small t/D ratio for a given material. 

Measurable deflections which are a small part of t force larger D. 
Fortunately, sensitivity is proportional to the fourth power of diam- 
eter. Conversely, attempts to reduce diameter are penalized by this 



ratio and the approximation. Timoshenko reports O .583 for ji = 0.25 
and 0.488 in another approximation. Our constant happens to be an 




D ~ 50 5 100 * 
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power© An attempt to achieve high sensitivity at given limited 
diaphragm diameter by reducing diaphragm thickness t (i.e.* small 
t/D) will ultimately face this limitation when the detectable de- 
flection is of the order of t© Another problem* described in the 
next section* is a less in sensitivity characteristic of thin dia- 
phragms© 



8o 3^ The Large Deflection of a Thin Diaphragm with 
Zero Initial Tension; the Solution of Hencky© 

The following comparisons of thin diaphragms* or membranes* 
with circular plates having bending stiffness will use the center 
deflection as a basic measure© It is well to observe that this 
may be an incomplete criterion* depending on the means employed to 
detect diaphragm distortion# Center deflection is sufficient if 
the detection system depends only on this* or if the deflection 
form of the cases compared is similar# The deflection form of stiff 
plates and membranes are not geometrically similar© Therefore the 
volume sensitivity* the volume swept by the deflecting diaphragm 
per unit pressure, is different in each case© A practical case in 
which this would be significant is when the change in capacitance 
is used to detect deflection© A membrane* sweeping out a larger 
volume than a stiff plate at the same center deflection* produces 
significantly greater unit change in parallel-plate capacity (see 
for example Lilly* Legallis* and Cherry [19^7])® Since the primary 
index is center deflection* attention will be given only this item. 

One may suspect that a sufficiently thin diaphragm would have 
no effective bending stiffness* so that deflections would develop 
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only membrane tension stresses. Kirchhoff theory and the solution 
of Way do not apply. H. Hencky [1915] has achieved an exact solu- 
tion by power series methods for the equations of equilibrium of a 
flat disc uniformly loaded on one side, subject to the assumptions 
t^ << t and no initial tension. 



A commonly used plate modulus defined as 




Et' 



12(1 - ji ) 



is in effect reduced to negligible value by Hencky *s assumption. 



Hencky’s result for deflection at the center of this disc is 



| = 0.662 




[The constant 0.662 is given by Hencky for Poisson ratio jx = 0.3; 
other authors have values 5 percent or less]. 



This result is essentially predicted by the approximate solu- 
tions for the large-deflection case discussed above with the added 
assumption 

t 3 6(i> o )2 t 

(i) 3 < < 0.545 C-n*- • 



so that from those approximate solutions 



tfF = 6(p o ) 
D = P o 




0 . 5^5 (• 



6(p 



-) (£) 



or in equivalent form 
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= o.68 




Hencky’s solution has the elegance of exactness; this coincidence 
is given as a suggestion of reasonable continuity between solutions 
recorded here. 



Hencky’s solution is also shown graphically with the pre- 
ceding solutions in Figure 35. 

The instantaneous sensitivity is conveniently found from the 
Hencky solution by logarithmic differentiation, 



so that 



d6 1 dp 
6 ~ 3 P 



d6 _ 1 6 
dp 3 P 



1 

3 



o . 662 




, P 



The initial sensitivity is indicated to be indefinitely large, 
which is apparent from the slope of the deflection curve. Unfortu- 
nately the very high initial sensitivities appear to be not realizable. 
A physical interpretation of the infinite initial slope of the deflec- 
tion curve may be that the no-load rest position of a thin diaphragm 
without initial tension is unstable. This is possibly a part of the 
phenomena called "oil-canning." Further, if one draws on the analo- 
gous result from the approximate solution, the necessary condition 

t 3 6(p o ) 2 t 

(— ) << 0.5^5 (—5 ) ( D ) 



is not fulfilled for practical (t/D) at sufficiently low 6(p q ). 
The assumption for this case implies large (D/t), small 
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t/D f so that very small pressures are equivalent to values of 

k k 

pD /Et producing operating points of 6, p away from the ini- 
tial region. A small hydrostatic preload does this for instance. 

For finite deflections the instantaneous sensitivity decreases rapid- 
ly to values less than the constant value of the ideal disc or the 
corresponding values of the large deflection (Way) solution. The 
theoretical average sensitivity of the Hencky-type diaphragms would 
appear to exceed that of the Way-type diaphragm. This is considered 
illusory because the instability of the no-load initial position of 
Hencky diaphragms would in practice mean an average sensitivity mea- 
sured from some initial non-zero deflection. 

This section has employed the Hencky solution to show that 
exploiting thin diaphragms without initial tension and in the limit 
having entirely negligible bending stiffness is not likely to lead 
to practical sensitivities equal to the theoretical sensitivity of 
the Kirchhoff diaphragm. A characteristic of this limitation of 
sensitivity is the assumption of no initial tension. We consider 
now whether there is benefit in adding initial tension. 

8.35 The Small Deflection of a Thin Membrane with Initial 
Tension. 

An exact solution to the equilibrium equations for small de- 
flections of a thin circular disc uniformly loaded on one side and 
subject to a large initial tension T is known. The formulation of 
this problem assumes that the bending stiffness of the plate 

_ -z p 

[F = Et /12(l»u )] is negligible. It is also assumed that the 
initial radial stress (initial tension) is sufficiently large so that 
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the pressure load is brought to equilibrium by opposing tension com- 
ponents due to the curvature of the membrane and not by additional 
stress® This solution is that of the membrane equation 



V 2 6(r) = ~ £ 

which yields for deflection at the center 




The deflection is linear in pressure; the pressure sensitivity is 
constant. 

In order to compare this with previous results, we make sub- 
stitutions, with s equal to the initial unit strain, 

o 



1 " % rrir 1 = % (1 + 



so that 



I = 255 ( F° ( 7 ) c 3(l^ u)c 'I 1 ] 






This is the same form as that for the ideal disc times the bracketed 
factor. That factor, for practical values of (t/D), reduces to 

_2f 

unity for initial strains e q of the order 10 or less. 

The pressure sensitivity of the membrane is constant 

b D 2 

« = F I5t ‘ 

In non-dimensional form, 



(5 F 3 
” = 25 ^ 
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If one compares this form with the result from Hencky’s solution, 
the factor 3(1 + ji)s Q /l6 and 0®5^5 (6/D) are seen to have equiv- 
alent roles in reducing the non-dimensional pressure sensitivity 
cfF/D * The apparent high sensitivities of the Hencky solution at 
small 6 are not realizable nor are those of the membrane solution 
for small initial strain since in either case the assumption inherent 
in the appropriate exact solution is not satisfied* 

We will now consider that the Hencky solution may represent 
one extreme condition of a diaphragm with no bending stiffness and 
the membrane solution another® The Hencky solution is for zero ini- 
tial strain® The membrane solution is for such large initial strain 
that radial stress may be considered constant. The region between 
should be considered since both Hencky and Membrane Solutions give 
the illusion of initial sensitivities greater than that of an ideal 
Kirchhoff diaphragm. 

8o 36 The Large Deflection of a Thin Diaphragm with 
Initial Tension® 

The membrane solution is not valid for small initial tension 
and Hencky’s solution includes no initial tension as a boundary con- 
dition® An exact solution of the equations of equilibrium for a 

3 

flat circular disc„ assumed sufficiently thin so that t ^ <( t t or 
in effect so that the bending stiffness is negligible, and subject 
to initial tension or initial strain^ is reported in Appendix B® 

This solution bridges the gap between that of Hencky and the mem- 
brane solution® 



As a result of this solution 
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t 



K 



\IZ z 

ViAs 



Here K for a given plate is a function of initial strain and the 
pressure (that is, the finite deflection) and is shown graphically 
in Figure 37. For zero initial strain 



K(0) = 0.662 ; 

that is, Hencky’s solution is reconstituted. K(e Q ) decreases in 
value for increasing g q . For large initial strain one has the 
asymptotic solution 



K(e ) ~ 
o 



1 (1 - p) H /e p 2 d‘ 



as s 00 
o 



so that 



P D 



Ee 



P D 

ITT 



as e cd . 
o 



l6t 



1 - 



The membrane solution is asymptotically recovered as it should 
be for large values of initial tension. 



The average sensitivity is 



^3_ ( kL d 4 

a ^ 2 ' Et l6 
P 



3 2 

For initial tension zero the term K /p tends asymptotically to 



y? 0.66 2^ 
2 ~ 2 
p T-0 P 



as p 0 



but this is not useful on the basis of previous remarks because of 
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Deflection coefficient k as a function of 

INITIAL STRAIN e 0 (per unit) AND LOAD* 




FIGURE 37. 



small residual bending stiffness or because of initial tension 
essential for stability of very thin plates. In case of initial 
tension T not zero, the constant tends asymptotically to 
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^ t ED 2 
p2 T^O t3 256 



recovering the membrane action. 



as p 0 , 



8.37 Composite Approximate Solution for the Finite 
Deflection of a Clamped Edge Diaphragm under 
Uniform Load. 

It remains to state a composite approximate solution for the 
finite deflection of a diaphragm taking into account bending stiff- 
ness, the influence of finite deflections and admitting initial ten- 
sion. 



Griffith [19273 reports carrying out an iterative procedure 
starting from the Kirchhoff (exact) solution for infinitesimal de- 
flections in the general equilibrium equations. His result in our 
algebraic form is 



tfF 

D 




(-) 

v d' 



,t, dm) (173 - 73P) ,5,‘ 
( D> *350 ^ ( D> 



(“) 

v d' 



5(1 + vO 



e 

o 



An argument will be made that this form, adjusted, is a correct one 
for an empirical correlation of the sensitivity of general diaphragms. 

The form reduces to the exact (Kirchhoff) solution in the ini- 
tial case with the deflection 6 and initial strain e q made zero. 

For initial strain 

o 



made zero and neglecting the cube of 
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t/D this i orm is comparable with the exact solution of Hencky and 
can be made coincident by small adjustment of the constant 
(1 + yu) (173 - 73 ju)/360. 

With only the initial strain term neglected this form is con- 
sistent with approximate solutions of the large deflection case 
which in turn are useful algebraic forms representing the exact so- 
lution of Way. The constant (1 + ji) (173 - 73yi)/360 may require 
adjustment but not to an extent practically conflicting with the 
preceding adjustment. 

For large values of initial strain and neglecting the dis- 

Z 

placement and (t/D)^ terms this form does not reduce to the 
exact membrane solution without a ten percent change in the constant 
5(1 + Ji)/24 to 45(1 + jx)/ 240 = 3(1 + ^i)/l6 and becomes exact with 
this change. 

Neglecting only the term (t/D representing bending stiff- 
ness, the sensitivity is controlled by finite displacement (Hencky) 
or initial tension (membrane action). 

The particular advantage of this empirical form is that it 
makes clear that the pressure sensitivity of a flat disc diaphragm 
will not be greater than the theoretical maximum of the Kirchhoff 
theory and displays that, for very thin diaphragms, the loss of 
sensitivity due to practically unavoidable values of 6/D or e o 
will be considerable. 

8.4 Dynamic Response of Flat Circular Diaphragms. 

8.4l Error Due to Dynamic Response Characteristics of 
the Diaphragm. 



The preceding investigation of the pressure sensitivity of 
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circular diaphragms was limited to static pressure loads and did not 
consider dynamic effectSo The dynamic response of the transducer 
should be considered if the pressure fluctuations include components 
at frequencies comparable with some fraction of the lowest signifi- 
cant frequency in the response of the diaphragm,* 



A convenient model for estimating dynamic error is to consider 
the pressure receiver to be a one-degree-of-f reedom linear mechanical 
system with a natural frequency equal to the lowest natural mode fre- 
quency of the diaphragm in free vibration® With this model the mag- 
nitude of dynamic error at pressure fluctuation frequencies u) small 

P 

in comparison with the diaphragm lowest-made frequency <o is pro- 

2 

portional to ( 00 ^/ 0 )^) and is minimum for a damping constant about 
70 percent of critical damping® 



For frequencies encountered in the experimental part of this 
work it is not difficult to obtain ratios (Up/(O n that are small© 

For other model sizes or more violent oscillations this may be more 
difficult® 

The high natural frequency of clamped stiff plates or stretched 
diaphragms is helpful in reducing this error® 



8okZ The Lowest Natural Mode Frequency of Circular Diaphragms. 

The dynamic response in vacuum of a circular flat diaphragm 
with clamped edges is studied in W® Mason [19^2]® The solution is 
applicable to diaphragms with given values of initial tension© The 
result is not conveniently expressed algebraically® A convenient 
graphical solution for the lowest mode frequency is given by 
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J. Patterson [1952] in the form CD n 2r 2 -ya 2 ) /t VEg as fun* 

2 2 2 

tion of 12e o r (1 )/t . This function is shown graphically in 



as func- 



Figure 58 « 

The least value of the lowest mode frequency for a given dia- 
phragm occurs for zero initial strain. For this condition the above 
form gives 



This corresponds to the clamped edge stiff plate equation for the 
lowest mode frequency found in handbooks as 



where V is the speed of sound in the diaphragm material* [The 
constants by various authors may vary less than five percent]. 

This form will probably be more convenient than the preceding be- 
cause of the practical difficulty of knowing E and, particularly, 

e . 

o 

A simplified form is also available for the case of large 
initial tension, corresponding to the membrane solution, 



The assumptions essential to the latter, or the complete 
solution first given, show that this value of 00^ must exceed that 
of the zero initial tension case* Sample calculations of the 
clamped-edge stiff-plate case show the conveniently high values of 




co n r/V = 3.0 t/r , 



uy/V = 2.4 VT/Et . 



0) 



n 



for small gauges if one takes the velocity of sound in thin 
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l2R 2 e 0 (l-p 2 )/ t 2 

GRAPHICAL SOLUTION FOR THE FUNDAMENTAL NATURAL 
FREQUENCY OF A CLAMPED, CIRCULAR, FLAT DIAPHRAGM 

( FOL LOW II N © PATTERSON [1932] ) . 



FIGURE 38 . 
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2 1/2 

plates (Eg/yo(l- ^ ) to be for example about 17,000 feet per 

second in steelo 

8.5 Pressure Gauges Mounted Internally. 

8.51 Advantage and Disadvantage 0 

The pressure sensitivity of a diaphragm-type pressure re- 
ceiver is a strong function of the diameter. The required pressure 
sensitivity is determined by the pressure fluctuation amplitude that 
is expected and the minimum deflection that is detectable. Further, 
the required diameter may exceed the maximum physical dimension al- 
lowable on the surface of the oscillating model. 

The apparent physical diameter on the body surface can be 
reduced by mounting the pressure gauge diaphragm internally. There 
is necessarily then an internal instrument cavity and connecting 
passage or tube. These introduce two effects which could require 
compensation in the measurement; they are (i) the effect of accel- 
eration of the enclosed fluid on the indicated pressure; (ii) at- 
tenuation of pressure fluctuations in the tube and cavity. 

8.52 The Effect of Acceleration of the Enclosed Fluid 
on the Indicated Pressure. 

The fluid enclosed in the instrument cavity and tube is forced 
to accelerate with the body. Therefore there is a pressure gradient 
in the enclosed fluid due to body acceleration. The top surface of a 
cavity is a suitable location for a horizontal diaphragm. 

As an example # consider an instrument cavity in the vertical 
flat side of a vertically oscillating wall. The fluid is assumed 
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An instrument cavitv in an 

OSCILLATING WALL. 



FIGURE 59. 
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inviscid, incompressible, and at rest. Because the wall and motion 
are vertical the fluid external of the instrument cavity remains un- 
disturbed. The external hydrostatic pressure is instantly communi- 
cated to the top of the instrument cavity by an infinitesimal port. 

The external hydrostatic gauge pressure at this port is p Q + y0g[6 y(t)], 
where 6 y(t) is the oscillatory motion and p Q is the hydrostatic 
pressure at the initial position. 

There is a pressure gradient in the fluid in the instrument 

2 2 

cavity due to body acceleration a = d y/dt • The pressure in the 
cavity is 

p c U) = yo(g - a)*l + P Q + fS C6y(t)3 , 

where tj is the vertical distance down from the cavity top to the 
horizontal plane at which the pressure is P c (n) • 

For a top port the indicated pressure on the top surface is 
a direct measure of the external pressure. 

8.53 The Dynamic Response of the Internally-Mounted Gauge. 

The question arises: what is the dynamic pressure sensitivity 
of the pressure receiver to pressure fluctuations at the entrance of 
the connecting tube? It is clear that the hydrostatic pressure sen- 
sitivity is not affected, but it is also well-known that dynamic ef- 
fects may appear. 

The analytic model considered here will exploit the conclusion 
found in Section 8 a that the lowest natural frequency of a flat cir- 
cular diaphragm will probably be high with respect to frequencies of 
interest in experiments visualized in this type work. We will now 
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assume that the lowest natural frequency of the diaphragm is also high 
with respect to the acoustic effects of the tube* cavity and resilient 
diaphragm* This assumption enables us to M uncouple M the dynamic re- 
sponse of the diaphragm* which accounts for diaphragm mass and elas- 
ticity* from the response of tube* cavity* and elastic (mass-less) 
diaphragm* Thus we expect dynamic behavior of the tube* cavity and 
elastic diaphragm to be of interest at frequencies low enough to con- 
sider the diaphragm to be stiffness-controlled* 

The problem is to consider the attenuation of oscillatory pres- 
sure fluctuations in a viscous fluid traveling through a connecting 
tube to an instrument cavity terminated with an elastic but mass-less 
diaphragm* This analysis is reported in Appendix C* Sinusoidal gauge- 
pressure fluctuations of frequency go are assumed at the entry port 
of the connecting tube of radius r* The pressure fluctuation at the 
diaphragm is assumed to be equal to the pressure at the cavity-end of 
the tube* The diaphragm elasticity and cavity volume provide a condi- 
tion on fluid continuity* The result is presented as a complex ratio 
of the indicated- to external-pressure fluctuation P^/P Q • The magni- 
tude of the ratio is the amplitude response normalized with respect to 
hydrostatic response* and the argument of the ratio is the phase lag 
of diaphragm sinusoidal response* 

The response P^/p Q is characterized by the value of the dimen- 

p 

sionless parameter 0 )a /V . This parameter also characterizes the ve- 
locity profile in the tube* (See L* Loitsyanskii [1957] * Section 86 * 

H. Schlichting [1955]? Chap* XI*) 

2 

For values of coa /-0 less than one, viscous forces are predominant 
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figure 40. 
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and inertial forces negligible. In this case 




1 + j(OT 



1 



where 



r = 8 /<d 2 a 2 
n 



C0 n = a characteristic frequency constant for a 

given fluid, tube, and cavity, and diaphragm 
elasticity. 



This result is also derived directly with a model assuming uniform 

Poiseuille flow throughout the tube and neglecting inertia effects. 

The relative response is down 29 percent at the cutoff frequency 

0) =■ 1/t and the phase error is 4-5 degrees lag. The cutoff fre- 

o 

quency is decreased with increased diaphragm elasticity (sensitivity). 



nificant. A convenient approximation is to present the system re- 
sponse in the form 



Here u is the normalized frequency ratio 0)/co 0 » where (D 0 is a 
characteristic frequency analogous to the resonant frequency of a 
moderately damped spring-mass system. For co = tU Q the indicated 



2 

For values of ooa /-J of about four inertial effects are sig- 




(1 - u 2 ) + ju 2 £ 



1 



pressure lags the driving source 90 degrees* The damping ratio is 

2 

For large values of 0)a /o the response ratio tends to 




1 - (co/(O n ) 2 



1 



This is useful in showing the undamped (non-viscous) nature of the 




Fractional response and phase 

LAG OP INDICATED PRESSURE, 
MODERATE DAMPING 

FIGURE 41. 
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acoustic response and also alerts one to the fact that strong acous- 
tic resonance is predicted which might damage a sensitive diaphragm.) 

The most general expression is 

£i _ aL 

p aL coshaL+ bL sinh aL 
o 

where 

aL = j 2 k(L/A) Vl/(1 - J) 
bL = - (0)/0) n ) 2 /(l - J) 

J = 2J 1 (ka)/(ka) J^Oca) 

(ka) 2 = jeoa 2 /^ . 

For tube lengths L small with respect to the acoustic wave- 
length in the tube A, this tends to 







1 

1 - J 



and approximations for 1/(1 ** J ) give the reductions cited above 
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APPENDIX A 

Mathematical Notes Referenced in the Text 



A-l Asymptotic Behavior of <p!J ( x ,y,h) for 

I x | -> cd, y = 0. 



A-2 



Asymptotic Behavior of cp (K,h) for 

s 



|x| a>, y = 0, 



A-3 Expanded Forms Useful in Numerical Calculation 

of Certain Definite Integrals. 

A-3.1 The Coefficients of q 2m^^ i- n 

Expansion Equations for the Method of Integration. 

A-3. 2 Two Special Forms. 

A-3. 3 Expressions in the Right-Hand Side of the 
Pg (6) Expansion Equations. 

A-3. 4 Alternate Forms for cp (6,©), Y„(&,©)» 

SSL S3. 

A-3. 5 Expressions in the Right-Hand Side of the 
q^mC & ) Expansion Equations. 
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APPENDIX A 



A-l Asymptotic Behavior of ^(x^.h) for |x| -»co, y = 0. 



The asymptotic behavior of the function 



QD 



f„ .. u _ x T /,-\--kh K smh ky - k cosh ky 

^p^Ax^h) = jt L(k)e “ : r—rr, “ coskxdk 

2m K coshkh- k sinh kh 

o 



for | x | ao , y = 0 as required in Section 3.32 may be determined 
by considering the complex contour integral 



x T r i \ K sinh ky - k cosh ky -kh ikx 
^ K cosh kh - k sinh kh e e 



dk = / F(k)e lkx dk , 



where the variable k is considered complex, temporarily. The 

integrand of cp" is an even function of x so that we may use 
cm 

| x | in the following. It is understood that h > Oo The contour 
of integration consists of R > Re k 0 , indented above the 
singularity defined by the positive real value of K q satisfying 



K cosh K h - K sinh K h = 0 
o o o 

2 

for given K = 0) /g and h > 0, the part of the imaginary axis 

R > Im k > 0 1 and the arc | k | = R, Re k > 0, Im k > 0. Here 
R > K q ; this is no restriction since we finally shall let R -*od. 
Integration around this contour yields 

/ F(k)e ikx dk + / F(k)e ikx dk + / F(k)e lkX dk + / F( k) e lkx dk = 0 , 

^ * 

with symbols indicating segments of the path. With the contour as 
defined , 
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ikx 



/ F(k)e dk -> 0 as R ao . 

Therefore after letting R ao and the radius of the indentation 0, 
we have 



/ F(k)is lkx dk = uiE t / F(k)e lk ^ x l dk 

f 



and, after taking the real part, 



00 I 

<p 2 m (x,y,h) = Re TtiR a + / [Re iF(ik) ] e" k 1 X 1 



dk 



where k is now a real variable and R is the residue evaluated 

s 

below. 



The residue R is 
s 



R = Lim (k-K )L e - kh k S - k C05k ^ e lk ' x l . 

S 0 K cosh kh - k smh kh 

o 



The limiting operation is executed after the substitution 



K = K sinh K h/cosh K h 
o o' o 



and manipulation of the terms 



(k - K ) 
o 



k cosh kh - k sinh kh 



(k-K ) cosh K h 

o o 

K sinh (k-K )h + (k-K ) sinh kh cosh K h 
o o o o 



With this substitution, the residue simplifies to 



-K h iK | x | 

R = L(K )e ° e ° K cosh K (h-y) . 

s o o o 



(k-K Q ) 

L:l “ k sinh (k-K )h + (k-K ) sinh kh cosh K h ’ 
k=K o o o o 

o 
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or, after evaluating the limit, 



Re niR = 
s 



- K „M K o )e 



"K> h - 

2K h sinh 2K h (h-y) sin K |x| 

O O 



We find 



■v , T , 2m 2m „ , r 2s 2s 

-K h K a , N a_ K a 

tc Trir 0 - -2 1 ^ v / 2n+l o 

- k ( 2m-l) 1 cosh K Q h (2e-l)| coshK Q h 



o o 



by using the definition of L(k) in Section 3.31(c) and substituting 

K = K q sinh K^h/cosh K Q h as well as evaluating the expression for 

k = K . 
o 



With this expression we may write 



cosh K (h-y) 

Re niR =• - E_ (K a,K h) ° ^ ■ sin K |x 

s 2m o ’ o cosh K h o' 



where 



E~ (K a,K h) 
2m o ’ o 



2n 



2K h + sinh 2 K h 
o o 



K 2m a 2m N 

O y / - vll C-iif- -L U 

_ (2m-l)( + n t 0 (2s-l)l 



^ 2s 2s- 
a_ , K a 
n 2n+l o 



The result for Re niR is identified as the asymptotic behavior 

s 

°f (p^ m (x,y,h) for |x| -» oo . On the surface y = 0, 



tp 2m (x ’°’ h) ~ - E 2m sin K o l X l aS 



X QD 



We note that E 2 m 0 as h oo . 
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A-2 Asymptotic Behavior of cp (K,h). 

s 



The asymptotic behavior of the function 



V. (K,h) = i g°3hk(h-y)coskx 

s K coshkh - k smh kh 

o 



for |x| -> qd as required in Section 3.34 is found by following the 
procedure of Appendix A-l. 



The contour integral to be considered is 



coshk(h-y) 

K cosh kh - k sinh kh 



The contour is taken to be the same as in Appendix A-l. After fol- 
lowing the procedure outlined there, one finds 

9 (K,h) = Re niR 1 + / [Re i F'(ik)] dk . 

s s 

o 

Thus the Re uiR* gives the first term in the asymptotic expression 

5 

for <p (K,h) as |x| -> oo ; that is, 
s 



2K h V sinh 2K h ~ C05h K o h C03h K ° (h - y) Sln K ° 1 x 1 • 
o o 

On the surface y = 0, as |x| oo , 

cp s (K,h) y=0 ~ E s (K o h) sinK o |x| 

where 

2 K (coshK Q h) 2 

E s (K o h) = 2K h + sinh 2K h * 
o o 

We note that E (K h) ^ n as h oo . 
s o 
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A-3 Expanded Forms Useful in Numerical Calculation of 
Certain Definite Integrals. 

A-3»l The Coefficients of p 2 m ( 6 ), q 2 m ( 6 ) i n Expansion 
Equations for the Method of Integration. 

Coefficients a_ , and b_ . are defined for use 
2m, k 2m, j 

in solving the expansion equations by the method of integration 
(Section 5 . 12 ). 



k/2 

/ f 2m (e) Sinks ae = a 2mA 



k odd : 



( 2 m+k-l ) 



2m( -1) 
4m^ - k 



2 



2 






6 

2m-l 



k even: k = 2m: - k/ 4 
k / 2m: 0 



k odd: 



< 



k even: 



k = m = 1: 0 

k = 1 / m: (-l) m k/ 4 
k = 2m-l / 1: + k/4 

otherwise: 0 

(, 2 m+k- 2 ) 

k 2m( 2m-2 ) ( -1) 

? 2 2 
k - 1 (2m-l) - k 



Y f 2 m (9) cos j©d© = b 2m< 

o 



j odd: 



2m 



,2 .2 T 2 m 

4m - 2 



= < 



2m-l 



2 m-l 



j even: -s — y 5 — p+(-l) 2 

2 - 1 1-3 2 J 

(zero is treated as an even number.) 



IT ( _1 ) m •< 



,i+i 1 

even: — j 



■ odd: r~ 

2 1+J 



m 1 



,i+2m-l 1 

“ — = even: ■= — z — r- 

2 2m-l- j 



,T + 2 m-l 1 

“ — - odd: := — = — r 

2 2 m-l+o 



1+2m 

iL -^ — even: 0 



+ -< 



i+ 2 m ... 4m _ 
2 ° dd: ta 2 - j 2 
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A-3.2 Two Special Forms. 



The following definitions are used for convenience 
of notation in later integral expansions. 

u/2 

I^Cl^n) = f sink© sin n© d© 

o 

0 if k or n are zero. 

= < + j- for n = k ^ 0 . 

for other k,n: 

(n+k) even: 0 

(n+k) odd: 



if k even: 



if k odd: 



( _ 1} (n + k + J)/2 

(_l)(n+k+ 1 )/2 . 



I 1 (n t k) 



tx/2 

/ cos n© sin k© d© 
o 



= 



0 for k = 0 
for other k,n: 



If (n+k) even, and if (n+k)/2 even, then 

if k even: 
k odd : 

if (n+k)/2 odd, then 
if k even; 
k odd : 



0 

l/(k-n) 



2k/(k 2 - n 2 ) 



l/(k+n) 



. 2 2 
k - n 



If (n+k) odd: 



148 



A-3.3 Expressions in the Right-Hand Side of the p 2ffl (&) 
Expansion Equations. 



*/2 


sin k© 


00 


/ *ca (6 > S) 


= n 


I (-1) 


o 


cos j© 


p=0 


m/2 


sin k© 





/ y c& (6,n/2) d© = n sin 6 

o cosj© I, (j # l) 



Pi I 1 (jfP) 
I 2 (k t l) 

1 < 



A-3.4 Alternate Forms for ip (6,©), y (6,©). 

S6L SSL 

The form of the quantities y (6,0) and cp (& t 0) 

sa sa 

used in the q 2m (&) expansion equation (Section 5»12) and in the 
calculation of M(6,Q) (Section 5.22) t respectively, are derived 
from the Cauchy Principal Value integral 



oo -ky +ikx 

* H rS" 

o 



£ -K(y-ix) Ei[K(y . ix)] 



oo n -in© 



= e -Ky [cosKx+ i sin Kx] [y + log r - i© + 7 - ■ , •] . 

- ° e ^ =1 n»n 



Here E^(-iz) designates the exponential integral of a complex 

argument [Tables of the Exponential Integral for Complex Arguments, 
National Bureau of Standards Applied Mathematics Series 51* May 1958]. 
Euler’s constant 0.57721 ... is designated y and © = arctan x/y. 

The real part of this expression represents a source in the free sur- 
face of fluid of infinite depth (Section 3.33) and the imaginary part 
represents the conjugate stream function. Contour integration leads 
to the real and imaginary parts given in Section 4.2, but the right- 
hand side of the above expression is one suitable to numerical 



^ o 



calculations. Evaluating this, 

^sa^ 6 ’®) = e 6c ° S ® [Q(9) cos ( 6 sin 9) + S(9) sin ( 6 sin 9) ] 
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- e CQ(9) sin ( 6 sin 9) - S(9) cos ( 6 sin 9) ] 



sa 



where 



CD , n 
0 



S(9) =9+7 ~~i — sin n© 

n=l n|n 



oo ft n 

Q(9) = y + log 6+7 — : — cos n9 

— e n _ 1 n'n 



It is noted that 



S(n/2) = n/2 - Si (6) 



Q(n/2) = Ci ( 6) , 



where Si (6) and Ci(6) are the sine and cosine integrals, 



A-3.5 Expressions in the Right-Hand Side of the q 2m (6) 



Expansion Equation. 



The following integrals appear in Section 5.12, 



it/ 2. cos 9 sin k9 

Y (6,9) sin 9 d9 - f e Q(9)sin ( 6 sin 9) sin 9 

sa o cos j9 



d9 - 



n/2 sink© 

- / e“ C0S S(9)cos ( 6 sin 9)sin 9 



d© 



oo . P ^2^’ P ^ oo n oo , ,p d' '* 

= (v ♦ iog c 6> K-d^ 1 §r ♦ Z Z c-d 1 ’* 1 

p=0 P * I-^jtp) n=l * p=0 P * I^jfP+n) 



cos j© 

I 2 (k,p+n) 



i 4 ,p 

* L c-D p+1 

p! IjU,p> • 



OD 

L 

p=0 



O 
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The quantities not previously defined in this appendix are 



IjCjiP) = \ CA( |p-j | ) + A ( | p+ j | ) ] , 

where 

( , 2 

(|*) for | x | =0 , 

|x| even, 
| x | odd; 



A ( | x | ) =J -^3- [-1 + (-1) X//2 ] for 



X 

V. 



and 



I^(k,p) = I* [B(|p+k|) + (signum k-p) B(|p-kj)], k 4 0 



- 0 for k = 0 



where 



0 for |x | = 0 



B( | x | ) ^ ^ (-l) (X+2)/2 for |x| even, 



~ (_l) (x_1) / 2 for |x| odd. 



x 



Tt/2 sin k0 u/2 sin k© 

Y oa (6,it/2)sin© d© = re cos 6 / sin© d© 

cos j© o cos j© 



sa 



co , -x u/2 sin k© 

* / £f— r / sine 

O 6 + X O COS j© 



dx d© 



I 2 (k,l) 

«= - C(rc/2 + Si 6)cos 6 - Ci 6 sin 6] 
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This form is found directly from the expanded form of y (6,u/2) 

SSL 

given in Section 5.11, and the identity 

CD , —X 

/ — f j dx = [n/2 - Si(6)]cos(6)+ Ci(6) sin (6) 

o 6 + x 

[Lamb, Hydrodynamics, Art 244(37)1 Here Si (6) and Ci(6) are 
the sine and cosine integrals and are available in tabulated or 
series form. Alternately, the same result for y (6, n/2) may be 
found in Appendix A-3*^ above where a series form suitable for com- 
putation is given. 
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APPENDIX B 



The Finite Deflection of Circular Plates with Initial Tension. 



B-l Introduction. 

B-2 Governing Differential Equations. 

B-3 Series Assumed for Solution. 

B-4 Evaluation of Coefficients b_ in Terms of b 

dn o 

and Definition of the Function c^ . 

B-5 Evaluation of Coefficients a- in Terms of b 

dn o 

and Definition of the Function c ^ . 

B-6 Evaluation of b Q by the Boundary Condition. 

B-7 Graphical Solution. 

B-8 Special Cases. 

a. Zero Initial Strain. (The Solution of Hencky). 

b. Lightly Loaded Plates.' 

c. Comparison with the Membrane Solution. 
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APPENDIX B 

The Finite Deflection of Circular Plates with Initial Tension. 

B-l Introduction. 

Consider a circular plate of homogeneous material obeying 
Hooke’s law which is initially flat until uniformly pressure-loaded 
on one side. The edge is clamped. The only widely-published solu- 
tion of the simplified equations of equilibrium including a given 
initial-tension boundary condition is the membrane solution. The 
membrane solution assumes vanishing bending stiffness, small de- 
flections, and that the initial tension is large compared to de- 
veloped material stresses. 

In order to maximize pressure sensitivity of a diaphragm- 
type pressure gauge, one is led to thin diaphragms of vanishing 
bending stiffness and to reduce the initial tension to the minimum. 
Application of the membrane solution is not wholly satisfactory 
when initial tension is small because material stresses become ap- 
preciable at moderate loads, and because deflections may not be in- 
finitesimal and because the solution gives no hint of departure from 
linearity. The exact solution of Hencky assumes zero initial tension 
and the infinite deflection derivative at zero load is not physically 
appropriate to stable diaphragms of finite pressure sensitivity. 

This appendix presents an exact solution to the equations of 
equilibrium for this type plate subject to arbitrary initial tension. 
For zero tension the solution of Hencky is recovered. For large ini- 
tial tension the membrane solution is approached asymptotically. For 
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lightly loaded plates the solution always tends to that of the mem- 
brane. The initial deflection derivative is limited by the initial 
tension, physically appropriate to the case of a lightly loaded, 
thin, stable diaphragm. 

The radial and tangential stresses, as well as deflection, 
are given in this solution. 

The membrane solution (large initial tension) is valid only 
for infinitesimal deflections. This solution is not so limited. 

The deflection pressure relation of the membrane solution is linear. 
This solution reveals the nonlinear solution, tending toward a 
Hencky-type solution as the developed stresses increase relative to 
initial stresses. 



B-2 Governing Differential Equations. 

The equations of equilibrium and stress compatibility for an 
initially flat circular plate with clamped edges were given by 
A. Foppl [1907] and are repeated by H. Hencky Cl915] and S. Way [1934] 
among others. 

We repeat the equations expressing horizontal equilibrium, 



& <V> 



■ 0 • 



and compatibility, 



E / d6 > 



dr 



<»r * °t ) * t <JF> ‘ 



- 0 , 



using <f and O', as radial and tangential stresses and 6 as de- 

r 



flection 
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The equation expressing vertical equilibrium, 



E 



1 - 



d /I d , d6 u pr . . db 

2 12 dr r dr (r dr'' ~ 2 + t0 r dr ’ 



is simplified in the case < < t, 



_ pr . . d6 

0 = “5r- + tcT — . 

2 r dr 



This is the case of vanishing bending stiffness. For extended argu- 

If If 

ment see, for example, Foppl or Hencky [ loc cit] or A. Foppl [1924]. 
We deal here only with the case of vanishing bending stiffness# 



B-3 Series Assumed for Solution. 



A power series form is assumed for the dimensionless radial 

stress ratio (o'/ E). The corresponding series for the tangential 
r 

stress ratio (tf^/E) is derived with aid of the equation expressing 
horizontal equilibrium. That differential equation is thereby satis- 
fied# 



= c. 



00 

1 • 

n=0 



2n 



<?> 



2n 



The stress is made non-dimensional by Young’s Modulus Ej the 

radius r, by the plate radius a. The dimensionless coefficients 

b_ are to be determined and c n is a dimensionless function of 
2n 1 

the uniform pressure load found later to be 

c i = b ( et ) • 

From horizontal equilibrium, the series form for is 
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<f. , <s 

Jl = r A. ( _£) 

E dr K E J 



+ 



6 

r 

E 



= c 



00 

lie 



(2n + 1) b 



2n 



<r> 



2n 



A series form is assumed for the dimensionless deflection 
(6/a) compatible with the condition 6=0 at r = a. 



6 

a " C 2 



CD 



n=0 



2n 



cd 2n' 

- £ a 2n <£> 

n=0 



The dimensionless coefficients a_ are to be determined 

kn 

and c^ is a dimensionless function of the uniform pressure load, 
later found to be 



c 



2 

2 




The center deflection at r = 0 is 

CD 

6(0) = (ac 2 ) J_ a 2n . 

n=0 



B-4 Evaluation of Coefficients b_ in Terms of b and 

d.n o 

Definition of the Function c^ . 

With aid of the equations of compatibility and vertical equi- 
librium we are able to express b 2n ^ erms after making 

a suitable choice of definition of c^ . 

The nonlinear equation of compatibility in the form 




d , T l \ JL_ / X UU' 

dr ^ E + E ' + 2 k E dr ; 



0 



with the equation of vertical equilibrium 
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°r d& _ 1_ ££_ 

E dr ‘ ‘ 2 Et 



takes the form 




+ 



1 

F 




o . 



The first term of this expression is rewritten as the product of 
two appropriate series, using for one factor 



r 





£[l 

Ln=C 



b (-) 

V n v a 



2n 



CD 

Z (2n+l) 
n=0 



b 2n< 



2n-i 

] 



so that 



= c. 




CD 

Z (2n+2)(2n) 
n=l 



b 2 „<r> 



2n-2 



* 



*>T [ l 



(2n+2) (2n)b 2n (“) 



2n-2- 



1 

F 






It is appropriate to define 

c i 3 = ( eT } 

in order to leave a non-dimensional expression which is used below 
to generate recursion formulae for b 2n • 

2 

r 00 2n_. r od 2n-2-. 

[ 2 b 2n<f> ] [ I < 2 ”* 2 ><2»> » 2n (|-) ] ♦ 8 . 0 . 

l-n=0 J L n=l - 1 

This series is expanded and coefficients of powers of r/a are set 
to zero individually, the first reads conveniently 

8 b-b +8=0 or b_b +1=0 . 

2 o 2 o 
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The succeeding equations for the next higher powers of r/a are 

3b, b 2 + 2b b 2 . 

4 o o 2 * 

6b 6 b Q 2 + 8 V>2 b 0 + b 2 3 = 0 * etc ' 



Solving these in succession, all succeeding coefficients b~ are 

2n 

expressible in terms of b . The first five are 



, to be determined from boundary condition 
b 2 = . i/b o 2 
= - 2/3b o 5 
b 6 = - 13/l8b o 8 
b g = - 17/lSb^ 1 
bi 0 - - J7/27b 0 11 ' . 



B-5 Evaluation of Coefficients a_ in Terms of b and 

2n o 

Definition of the Function c ^ . 

With use of the equation expressing vertical equilibrium we 

are able to find all coefficients a_ in terms of b after making 

2n o 

a suitable definition of the function . 



The form of the equation of equilibrium 



_J 1 d6 1 ££ 

E dr 2 Et 



with the series for cs^/E after using 



3F " - *2 Z «« 2 >»2n 

n=0 



CD 



2n+l 



oo 2n 

r ^ / _ \ ,r , 



- 2c ? - Z. (n+ 1 ) a, (-) 

2 a 2n a 

n=0 
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is 



ao 2n cd 2n , 

2 V2 ZW? Z n (a ' >1) a 2n<?> 

n=u n=0 



Substitute the previously selected definition of ^ and define 



c 



1 pa 1 
2 " £ Et c 1 




in order to leave a dimensionless expression which can be solved 
successively for ^ 



co 



I b 

n=0 



2n 




go 2n 

Z tn+l)a (f) 
n=0 



1 = 0 . 



When this series is expanded and the coefficients of powers of r/a 
set equal to zero individually, one has a series of equations 

0 = a b - 1 
o o 

0 = 2a_b + a b_ 

2 o o 2 

0 = 3a i( b Q + + a Q b ^ , etc. 

Substituting for b 2 n in terms of b Q and solving successive 
equations yield for the first six a 2 n 



a o = + 1/b o 

a 2 = + V2b o 4 

a l + = + 5/9b 0 ? 

a 6 = . 55/7?-b o 10 

a fi = + 105/90b 
o o 

a l0= + 205/108b o . 



l6o 



B-6 Evaluation of b Q by the Boundary Condition. 

The unit strain e in terms of the in-plane stresses is 

°t ’ = E£ * 

The edge r = a of the plate is rigidly clamped. If the plate 
is not initally strained, e = 0 at r = a and must remain zero be- 
cause of clamp action. If the plate is initally strained and then 
clamped at r = a, the strain e = at r = a and by clamp action 

must remain e Q . The boundary condition is 

cr - p* T - Ec 3 , 

r=a 

e = initial strain, 
o 

The series form for cf^/E and cf^/E at r = a are 

^t 3 = C 1 E l n (2n+1)b 2n 

r=a n=0 

V - £ b . 

r=a n=0 

The boundary condition in series form is 
ao 

I < 2n+1 " f )b 2n = e o /c l 
n=0 



or, expanded, 

b o + feS b 2 + b 4 + b 6 + * ‘ ‘ = e o /c l (1 -/ i) 



and in terms of b 



(3-/u) 1 2 (5-fO 1 13 (7-;*) 1 

' (1-p) b 2 ' 3 (1-/) b 5 • IF (l-p) ^F 



e o /c i (1 7 i) 



l6l 



For given material and load this equation can be solved for a 

numerical value of b and the problem is solved. 

o 



B-7 Graphical Solution. 

A convenient procedure for evaluating b Q from the boundary 

00 

condition equation and the deflection coefficient K = 21 a ? is a 

n=0 

graphical one. 

The left-hand side of the boundary condition equation, includ- 
ing terms up to the n-power in b Q , is designated f n (b Q ) so that 
with sufficiently large n and b Q > 1 it is nearly true that 

W = £ </ c l (1 • 

Values of b Q were assumed in the range 1 < b Q < 10 and n was 
taken as 17. The numerical values of f^Cb^ are Pitted i n the 
figure to enable graphical solution of the inverse problem: given an 

initial strain e q the quantity e Q /c^(l - yu) is computed and the 
corresponding b Q may be determined graphically. 

The center deflection is 



CD 



6(0) = (ac 2 ) L a 2n = K T Et - 



n=0 




where we have introduced K, defined by 



CD 



K - £. a 2n = f - + rr * 

n=0 o 2b 

o 



5 + -^r + 



9b 7 ?2b 10 

o o 



70 205 

-1 5 + T C + • • • • 

60b 15 108b 15 

o o 



This summation has been computed for the assumed values of b Q , 1< b^ ^ 10, 
using these six terms, and the numerical value of K is presented 
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ThF_ FUNCTIONS ^ 7 (b 0 ) AND K(b 0 ) 



FIGUPE. 42 
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graphically. As a matter of convenience, the two figures 

f 17 ( V vs b o 



I a 



n=0 



2 n 



vs b 

o 



are graphically combined to eliminate b and present 

o 

K vs £ o /c 1 (l-yu) . 

With this graph, one enters with known load and initial strain to 
find a value for K and with this the center deflection is 




B-8 Special Cases. 

a. Zero Initial Strain. (The Solution of Hencky). 

If the initial strain is zero, this solution be- 
comes that of Hencky. The series expression f(b Q ) is equal to zero. 
The corresponding value of K is evaluated for the root b Q of 
f(b Q ) = 0. Our computed value using f^y(b Q ) differs by less than 
one percent from Hencky ’s value K = 0.662 and the difference is no 
doubt entirely due to numerical procedures. 

b. Lightly Loaded Plates. 

This special case is appropriate for considering 
the initial pressure sensitivity -and for comparison with the zero- 
initial tension case which predicts an infinite initial-deflection 
derivative suggestive of an unstable initial position. 



As the pressure load approaches zero, c^ ap- 
proaches zero, and f(b Q ), proportional to l/c^, is increasingly 
large. The function f(b Q ) approaches its asymptote, b Q : 



f(b )] 
o 



b^ large 



= t>-i£ 

O 1-JJl 



b * 

= o 



For large values of b Q , the value of K is 



K 



00 



Z a 
n=0 



2n 







~ 1 



" b 



o 



- 1 
= fTTT * 

o 



One may substitute for f(b Q ) its equivalent in 
terms of the initial strain, g q /c^( 1 -p), but it is more revealing 
to make another substitution as follows. 



Under the assumed conditions, b Q is large and 

therefore all b^ n for n > 0 are negligible. Therefore, at the 

edge r = a, o' = o'. and the edge boundary condition may be written 
r t 

*t ’ = Cr(1 "/ l) = Ee 0 ] r=a * 

The initial tension T = eft, so that one may use for f(b Q ) an ex- 
pression in terms of T: f(b Q ) = e^/c^Cl-p) = T/c^tE. With this the 

center deflection coefficient K is c^E/T and the center deflec- 
tion reduces to 

6 ( 0 ) = ^ 



which is exactly the membrane solution 
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Thus the initial deflection of a lightly loaded plate is con- 
trolled by the initial tension. It is not considered appropriate to 
assume zero initial tension for a stable diaphragm of vanishing bend- 
ing stiffness; the initial deflection derivative is bounded and equiv- 

2 

alent to the membrane result a /4 t * 

Co Comparison with the Membrane Solution. 

The membrane solution may be recovered here by 
observing that a particular solution to the equation expressing hori- 
zontal equilibrium is 

= cr. = cT = constant = T/t . 
r t 

The equation of compatibility will be satisfied only for in- 

2 

finitesimal deflection in this case, so that (db/dr) is negligible. 
The vertical equilibrium equation is equivalent to the membrane equa- 
tion , 

y 2 6 = - P /T . 

In terms of our series solution, b^ n are zero for n > 0 and the 
result is given in section b above. 

The complete series solution given here is not limited to in- 
finitesimal deflection and shows that as the pressure load (deflec- 
tion) increases. c„ increases, f(b ) decreases and the deflection 
T 1 o 

tends to that of Hencky [f(b ) = 0] as the developed stresses be- 

o 

come increasingly large compared with the initial stress. 
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APPENDIX C 

The Attenuation of Oscillatory Pressure Fluctuations Travelling 
in a Viscous Fluid Through a Connecting Tube to an Instrument 
Cavity with Resilient Walls. 

C-l Description. 

The required diameter of a diaphragm-type pressure gauge may 
be controlled by the sensitivity desired. This diameter may exceed 
the allowable diameter on the surface where pressures are to be mea- 
sured. The apparent diameter can be reduced by communicating the 
pressure through a port of allowable diameter to a cavity whose one 
wall is the required diaphragm. 

The apparent diaphragm diameter is reduced, but there may be 
a penalty in degradation of the dynamic response of the system tube- 
cavity-diaphragm compared to that of the diaphragm alone. 

This appendix reports an analysis of the attenuation of oscilr- 
latory pressure fluctuations travelling in a viscous fluid from an 
entry port through a connecting tube to an instrument cavity. The 
instrument cavity is expected to have a resilient wall, the diaphragm 
of a pressure gauge in our application. Related work is known in the 
literature but is usually concerned with long tube lengths, compres- 
sible gases, or does not include the effect of a resilient wall. 

The present analysis concentrates on the arrangements expected 
in measurements of small pressure fluctuations on the surface of 
small bodies oscillating in a fluid. The slight compressibility of 
a fluid and the high flexibility of a sensitive diaphragm contrast 
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with the more usual situation of gas compressibility and relative 
cavity rigidity reported in acoustics. Diaphragm elasticity plays 
an important role. It may in fact be the limiting factor in dynamic 
measurements of this type. 

C-2 The Momentum (Navier-Stokes) Equation Yields Particle 
Displacement and Volume Rate of Flow. 

The momentum equation for a fluid element in the tube is the 
Navier-Stokes equation for the axial velocity component w . We as- 
sume radial symmetry. 
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d w 



dw 



a . + w a 
d t dz 



A 2 
. d w 



dr 



1 — 

r dr 



a 2 
d w, 



. s. *£. + ^ (i-i + i — + ) 

yO dz V ° - a - 



dz 



Assume that dw/dz is negligible compared to dw/dt and dw/dr. 
Assume also that the particle displacement is harmonic in time* Then 



6(r,t) = 6(r)e*^ U)t , 
w(r)e ^t = |S.. Jtt6(r)e j«t t 

. + V . i 2* ) = £ 

dt 2 r dr' <o dz * 

dr r 

d 2 6 1 dfi 0)6 8_ d£_ 

dl> 2 + r dr + jV jo>ji dz 

A particular solution of this equation, assuming p not a function 
of r, is 



6 = (ffOs/CD 2 ^ . 

A homogeneous solution of this equation, finite at r = 0, is 
6 = J Q (kr) 

k 2 = - jo )/-J . 

Combining these solutions with a constant c, 

2 

6 = (dp/dz)g/u> yO + c J Q (kr) . 



The constant is evaluated by the no-slip condition that 6=0 for 
r = a . 



(£L )k 

K dz JS 



2 

0) p 



J o< to > 1 

J 0 (ka) 



6 



1 
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It is noted that the radial derivative of particle displace- 
ment f d6/dr t vanishes at the centerline r = 0 as it should. 



The volume rate of flow through a given section is found by 
integrating the area-velocity product over the tube radius at that 
section. 



dV = 2*r dr = [ r dr . 

/MO |_ J Q (ka) 



3t 



] 



V = 



a 

/ dV 
0 



,12ng(3p/3z) 



yOGO 
2 

,1na g(3p/3z) 
j ogo 



2 k 2 J (ka) 0 
o 



ka 

/ kr J Q (kr) d(kr) 



2J x (ka) -J 
kaJ Q (ka) 



1 - 



For simplicity we use notation 



J = 2J 1 (ka)/kaJ Q (ka) 



V = - p( dp/dz ) 



>G0 



P = na g(l-J)/jyOC 
It may be observed that for small (ka) 



1 - J (ka) /8 . 



Substituting this in the expression for the volume rate of flow, 

V = - g( 3p/3z)/8 , 

the known volume— rate of Poiseuille pipe flow in which inertia forces 
are negligible. That is, for go /-J vanishingly small this solution re- 
covers the Poiseuille solution. 
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C-3 The Equation of Continuity. 



The equation of continuity of fluid in a thin section of a tube 
of area A and flowing with radial symmetry is 



We will substitute the previous result 
V = wA = - p( dp/dz) 

and change variables from the actual local pressure p to the frac- 
tional overpressure tp . The overpressure cp, or pressure signal, is 
normalized with respect to the amplitude p Q of the sinusoidal pres- 
sure assumed to exist at the tube entry. The forcing pressure p Q 
is assumed to be independent of r . 



C-k The Equation of State and Effect of Tube Wall Elasticity. 

Introduce now the equation of state of the fluid in the tube 
in the form of a relation yo(p) . For the small pressure fluctuations 
expected in this work the density-pressure relation may be expressed 




T = (p - P 0 )/P Q 



dz *0 dz 




by means of a constant bulk modulus of compression of the liquid: 




Analysis from this point concerns nearly incompressible liquids. For 



gases, a similar path would be followed but different equations of 
state are necessary (A. Iberall [1950]). 
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The elasticity of the tube walls may not be negligible in com- 
parison with the slight compressibility of the liquid. In this case 
the cross-section area A is a function of pressure and might be 
written as approximately 

P - P„ 



A = A (1 + f 
o 



E 



-) 



where E = Young* s Modulus of wall material, 

f = a geometric function of wall thickness and nominal 
form. For a thin-wall circular tube f = a/t , 

A — cross-section area at p = p . 



[For a table of f for ordinary shapes see I.P. Ginsberg, Prikladnaya 
Gidrogazodinamika, Leningrad University, 1958, p. 106.] 

Combining these results to define an apparent bulk modulus in the 
tube, we find 



P k =/> < 1 + 



l - V . „ f(p - V 

o' B, ^ A o ^ + 



Pc k_a ♦ 



o o 



1 

P - P c 
B 



L ) , 



) 



B B 1 + E * 



Whether this apparent modulus B differs from is a measure of 

the effect of tube wall elasticity. Using the above expression for 
(^oA), we obtain 

d/>A _ P o^o^o dip 

at ~ b at * 
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With this result, the equation of continuity is 



p A 

IOO 



&£. _ ,os *L± = 

at r p 2 

dz 



= 0 



Here we have required that the product yO(3 be independent of z . 
It is noted this product is independent of density, so that we have 
not assumed density independent of z but we have assumed that the 
kinematic viscosity is sensibly independent of z . 



C-5 Boundary Conditions and the Effect of Diaphragm Elasticity* 
The boundary conditions are these* 

(i) The pressure fluctuation at the tube entry is specified. 



cp = <p Q e 



jcotj 

z=0 



(ii) The mass rate of discharge from the tube equals the 

accession of liquid in the cavity. This augment to the 
contained fluid is accommodated by elasticity of the dia- 
phragm and compressibility. 

The second boundary condition is expressed analytically as follows. 



The fluid mass in the cavity is jo 



V 



ment is 



dv . 

o .. do 

/°o dt + o dt 



From the previous relation yo(cp) we derive 



so that the mass aug- 



dp _ Po_ d£_ 
dt " B 1 dt * 

To determine the rate of change of instrument cavity volume we in- 
troduce the volume sensitivity, s, of the diaphragm, 
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dV 

o 

j — S • 
dp 

The volume sensitivity is geometrically related to the pressure sen- 
sitivity cT of the diaphragm* a measure of the change in center de- 
flection per unit change in pressure. For a membrane diaphragm de- 
flecting to a paraboloid of revolution, 

s - o' nr^/2 , 

With the volume sensitivity defined as above, we find 



dV ^ 

o dp 

dt S dt 



The mass influx into the cavity is 



£° * . s is. = (V „ SB. ) 

_ dt / o dt B, o 1 dt 



The expression (V q + sB^) represents an equivalent cavity volume 

V corrected for the resilient diaphragm. The term y° 0 /B^ is 
2 

equal to g/c , the squared velocity of acoustic propagation in 
the volume. With these notations the mass influx 



*2 ft * 



is equated to the mass flux at the end of the tube, 

/° v] = - ^ fJr ] = - ff ] * 

z=L z=L z=L 



The second boundary condition is 



0] = . v = - , 
0 Zj T q 2 

Z = L yOpC 
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where 

b = jo)g V^oPc 2 . 

C-6 The Solution of the Governing Equation. 

The differential equation from Section 3, after combining the 
equations expressing momentum conservation, continuity, and the ap- 
parent equation of state, is 
,2 A P 

0 C£ _ o ' O Q(p 

3z 2 = />?,* at • 

The solution is 

, az -az x jcot 

cp = (c-^e + c 2 e ) e J 

where 2 

a = . 

The desired form of the result is cp/ c p which is 

c c 

t-i . r az ] / t-i + i: . 

C 2 C 2 

The ratio c ]_/ c 2 I s eva -l ua ted by means of the boundary condition of 

continuity at the cavity, cp + bcp = 0 at z = L, so that 

z 

c n . v -aL 

1 _ (a-b)e 

c, “ , , \ +aL 

2 (a+b)e 

Finally, evaluating the pressure-signal ratio at z = L, 

(aL) 

cp ” (aL) cosh (aL) + (bL) sinh (aL) 
o 

The ratio cp T /cp is the pressure fluctuation in the cavity per unit 
L o 

fluctuation at the tube entry. This ratio is called Pj_/P D i n the 
text. The products aL and bL are as follows. 
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aL 



. Ml 1 
c t VlTT * 



Here , the apparent sound velocity in the tube, is yO Q /Bg, 

determined by the apparent modulus B rather than the fluid mod- 
ulus alone as in the cavity. The ratio oo/c ^ is the fluctuation 
wavelength in the tube, i. e . 



aL = j 2it ^ 



* 

Similarly, for bL we have 



bL = - (o' 



2 . 



2 VL 



2 a 1 - J * 
c A 



The quantity c A/VL is the square of a characteristic frequency, 



wj; = c 2 a/vl , 



so that 



bL = 



O) 



< 0 _ 



2 1 - J 



It should be noted that aL and bL are, in general, not 
real variables, as J is the ratio of a complex quantity and Bessel 
functions of a complex argument. 

The natural frequency co n is lowered by increasing diaphragm 
sensitivity through the effective volume V . 



C-7 The Short Tube Approximation. 

If the tube length L is a small part of the pressure fluc- 
tuation wavelength in the tube, the product aL is small and 
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P i aL 

P Q (aL) cosh (aL) + (bL) sinh (aL) 



- 1 



1 + bL 



1 - C-fiL) r-i 
0 ) 1 - 
n 



This is the form discussed in the 



J 

text. 



C-8 



2 

Approximation Valid for Small Values of (ua /-J . 
The value of 



2 J 1 (ka) 

^ ~ ka J (ka) 
o 

for small values of the argument 
(ka) 2 = - jcoa 2 /-J 



is approximately 

J 



1 + 



(ka) 2 

~~8 



so that 




The response ratio. for this case is 




1 

1 + j(l)T ’ 



x = 8-«yu) n 2 a 2 . 



This result is identical with the one obtained by assuming 
uniform Poiseuille flow and constant pressure gradient. If the 
Poiseuille flow mass rate, 
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Q V v --*%■* (f|) . 

with constant pressure gradient, 



d£_ P i - P Q 
dz L 



is equated to the instrument-cavity mass rate 

gV dp i 

2 dt ’ 
c 



the governing ordinary differential equation is 



dp i 

T —* p i = p o (t) • 



Assuming periodic flow yields the result given above. 

The time constant r is increased by increasing diaphragm 
sensitivity or decreasing the entry port diameter, for 

T = 8VLV/rca^c 2 . 



C-9 Approximation Valid for Moderate Damping Values. 

The identity 

1 J o (ka) 

1 - J J 2 (ka!) 



is, for moderate values of the argument, approximately 

k . 8 

•2 ” 3 p • 

^ ( cua /-J ) 

The corresponding response ratio is 

£i = i 

P o (1-U 2 ) + ju 2 t, 
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where we have used the notation 



u 





This form is motivated by analogy to that for the response of a 
moderately damped spring-mass system. This system responds with 
90 degree phase lag at a dimensionless frequency u = 1.0 with an 
amplitude ratio 1/2 £, . 

This form is also analogous to the classical Helmholtz 
Resonator results, but the present more-detailed analysis predicts 
a lower resonant frequency of 90 degree phase lag and generally 
greater phase lag at frequencies near resonance. 

To clarify the differences in the results, compare the above 
with the following. We recall from section 8 the ordinary differ- 
ential equation 



Multiplying by the tube area makes more evident that this is a 
force balance of pressure force and Poiseuille tube-wall friction 
considered to be constant along the tube length: 



tube-wall friction per unit area. Add to this force balance an 




na g 

where (4VQ/na^g) is recognized as the constant Poiseuille-f low 




inertia term 



and the force balance may be written 



d 2 p. 2. dp. 2 . 2 

l c A c A _ c A 

dt 2 LV dt LV Pl LV P ° 



The steady-state solution, assuming periodic pressure fluc- 
tuations of frequency 0 ), is 



!.(“.) + J2 

"n % 



where 



0) 



n 



o 2 a 

2 c A 
LV 



as before and 



^o = t/2 = ^/ a \ • 



This is the result given by Rayleigh and generally observed 
to predict resonant frequencies that are higher than measured values 
(M. Colby [19393 p. 220} J. Patterson [1952] p. 9). 

The more detailed analysis, admitting longitudinal variation 
of pressure gradient and a more general radial velocity profile, 
predicts acoustic resonance at a frequency about 14 percent lower. 
This further lowers the expected range of frequencies free of acous- 
tic effects. It is observed that in either case the effective in- 
strument cavity volume is used in calculating co n . 
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C-10 Approximation Valid for Cases with Small Damping. 

The identity 

J (ka) 

1 - J J 2 (ka) 

tends to unity for large values of the argument. The response 
ratio is therefore 




The value of this approximation lies in the fact that it emphasizes 
the tendency toward undamped acoustic resonance. 
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